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Redresseurs synchrones thermiques 


Par JEAN-PIERRE Bore et Piet Cornaz, Lausanne?) 


1. Introduction 


Dans beaucoup de domaines de la physique et de 1’électrométrie, l’augmen- 
tation de la sensibilité des appareils de mesure est limitée par le phénoméne de 
bruit de fond. Celui-ci est généralement di au caractére discontinu de I’élec- 
tricité (tubes, photocellules, semiconducteurs etc.) et a l’agitation thermique 
des électrons des conducteurs. Ces perturbations sont inévitables et il s’agit de 
limiter par des appareils appropriés leur influence sur le résultat de mesure, 
c.-a-d. d’augmenter le rapport signal sur bruit. 

Dans tous les cas ot la tension 4 mesurer est sinusoidale pure ou s’il suffit 
de mesurer une composante déterminée de son spectre de Fourier, la détection 
-synchrone, appelée aussi détection de phase ou détection lock-in est un moyen 
‘simple et efficace d’augmenter le rapport signal sur bruit [1]?). C’est un des 
“meilleurs systémes, en particulier si ce rapport est inférieur 4 1 A l’entrée [2]. 
On utilise le fait que le signal est périodique de période connue, qu’il a donc 
un signe déterminé pendant des demipériodes connues a l’avance, tandis que 
le bruit, représenté par une fonction aléatoire est tantét positif, tantot négatif 
pendant ces mémes intervalles de temps. Pour tirer parti de cette différence 
-fondamentale, on module la tension d’entrée (signal + bruit) par une tension 
de référence périodique d’amplitude constante de méme fréquence que le signal | 
et de phase fixe. Comme nous le verrons plus loin (§ 2) le produit contient alors 
‘une composante continue proportionnelle au signal et des composantes alter- 
natives dues au bruit. Un filtre intégrateur permet alors d’éliminer trés effica- 
‘cement ces termes alternatifs et de ne garder que le terme constant contenant 

Yinformation du signal. 

_ Lamultiplication peut se faire en principe par tout élément non linéaire; on 
utilise couramment des tubes pentodes, triodes, diodes ou des diodes a cristal [3]. 
Le but de ce travail est d’étudier quelques types peu connus et particuliérement 
‘simples de redresseurs synchrones utilisant l’effet quadratique d’échauffement 
da au courant. Ce principe a été proposé et utilisé par l’un de nous (J. P. BoREL, 
en collaboration avec Cr. MANus) sous la forme présentée au § 3 dans un 


as 
__}) Dans le cadre du Subside du Fonds National N° 1447. Ecole Polytechnique de I’ Université 


ode Lausanne, Laboratoire de Physique Technique, directeur: Prof. R. MERCIER. 
S 2) Les chiffres entre crochets renvoient 4 la Bibliographie, page 101. 
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spectrographe de résonance électronique [4]. La variante du §4 en est une 
amélioration mise au point dans ce laboratoire au cours de nos travaux en 
résonance nucléaire. Elle a une sensibilité meilleure, tout en étant plus robuste ~ 
et plus économique. 

Remarquons qu’un principe un peu semblable est connu dans le domaine ~ 
de I’électrotechnique pour réaliser des Wattmeétres [5]. Effectivement, tous les — 
schémas de lock-in qui vont suivre peuvent étre facilement adaptés au probléme 
de la mesure d’une puissance ou d’un cos ¢. 


2. Principe général 


Le principe des redresseurs synchrones thermiques est le suivant: 2 éléments ; 
dont une caractéristique électrique est fonction de la température (thermo-_ 
couples, résistances métalliques, thermistances) et deux résistances sont montés 
en pont (figure 1, cas de thermocouples). Celui-ci est alimenté sur deux points - 


stig 


Figure 1 

Principe du redresseur synchrone thermique avec thermocouples. : 
4 
opposés par le courant d’entrée, mélange d’un courant «signal» a ae 
t, = I, sinw #, et d’un courant perturbatewr, par exemple un courant de bruit_ 
que nous teprésentons par ip = J’ Ip; sinw,; ¢. Aux deux points opposés on 
applique un courant de référence i, = I, sin(w t + p) beaucoup plus grand qu 
la valeur efficace des deux autres. 

A cause de la cohérence de phase entre le courant signal et le courant d 


référence la puissance instantanée P, dissipée dans la résistance 4 sera diffé- 


: 
. 


he 
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rente de celle P_ dissipée dans 7: 


0 fa 2 


: Boks ase, 
Pi =R=R (= sing t + 2) ae sinw, ¢ —- oF 


sin (w ¢ + )| 


qT? P “ + i ; . . 
=R [= sin2@ t + Dap - ar ae sinw, t sine; ¢ 
4 ] 


rP nats a) 
+ GE sin? (wt + ) + Day ae ®t sinw tsinw,; i 


v 


Med é ets : 
Sg 3 Sinw tsin(w t+ y) + oO eee sin (w ¢ + ¢) sine; i). | 


Le premier crochet représente une puissance P,(¢) fonction du temps qui est 
la méme pour les 2 résistances 7 et 4. Les thermocouples ne mesurant que les 
différences de température (montage en opposition) ne réagissent pas a ce 
terme. I] ne sera déterminant que pour Ja température moyenne. En général 
on prend I, $I, et I,;; il suffit alors de garder les termes en J? soit: 


JT? — [2 
Pye R [1—cos2(wit+9)]; RhaeR, (2) 
ot P, est la valeur moyenne (au point de repos). 

Le deuxiéme crochet dans l’expression (1) donne l’écart de puissance par 

rapport 4 P, qui produira une différence de température détectée par les 
: thermocouples. Il peut encore s’écrire en posant Jw; = w; — @: 
4 I, q, 
SPH R {Hs [cosp — cos (2 wt + q@)] 
3 (3) 
+ »' aso [cos (Aw; t — p) — cos ((2@ + Aa,) t + e)}. 


i 


Cet écart de puissance devient maximum pour y = 0 ou z (c.-a-d. signal et 
_ référence en phase ou contre-phase). Pour y = 0 on obtient 


1 


ay 


fet 


cos2mt 


tek leis 
Pa = LR et + Dy coset — 


he 


nis = Fr FBé cos (2 @ + Aa,) i\ ae 
7 ‘i a 


bi la 


Cette expression contient un terme constant proportionnel au signal, des 
termes lentement variables proportionnels aux J; et des termes de haute 


fréquence. 
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L’inertie thermique du filament joue le réle d’intégrateur. Sa température © 
suit le premier terme qui est proportionnel au signal et avec une certaine 
atténuation le second terme qui varie lentement ; les termes de haute fréquence 
ne se font pas sentir. Cette atténuation est d’autant plus forte que Aw; est plus 
grand ce qui ne laisse passer que les perturbations de frequence trés voisine de 
celle du signal. 

Pour analyser de plus prés cette atténuation on peut faire les remarques 
suivantes: 

Soit c la capacité thermique de l’élément et soit AP =a AT la variation 
de puissance nécessaire pour provoquer un écart stationnaire de température ~ 
AT par rapport au point de repos. L’équation différentielle pour l’écart AT en — 
régime non stationnaire est, en négligeant une résistance thermique éventuelle 
entre fil chauffant et élément sensible: 


Oe 4+aAT = AP(). (4) 

Dans les thermocouples travaillant dans leur région purement quadratique « 
est une constante indépendante du point de repos. Pour les lampes ou thermis- 
tances sous vide la fonction T = /(P) n’est pas linéaire. Conduction et rayon- — 
nement y interviennent simultanément et la relation est assez bien approchée ~ 
par une expression du type 

Ped(T=T)+¢e(14— 1); (5) 

ou. T, est la température ambiante. Le coefficient « qui mesure la pente de 
cette courbe au point de repos dépend alors de celui-ci. 

Pour AP(t) =a+ )D/ 6; cos(w;t + y,) Véquation (4) a la solution perma- 
nente: 


ri, a b; ' ; c 
A Doe of 2 re ae cos(w;t+9;—6;); avec tgd,; = 3%: ©) 

On voit que les termes de fréquence w sont atténués comme par un filtré RC 
de constante de temps t = c/a. Cette constante est de l’ordre de 1 s dans les cas" 
usuels. En prenant pour A P(t) l’expression P, de (3) et en négligeant d’emblée 


Ls 


les termes en w et 2 @ on obtient pour l’écart instantané de la température: 


] 
1 EOL I, Ip; il ; 

AT, =+ R [tt + te ee cos (Aw, ¢ — | ' 

ae a“ 4 2 4 y1 eee Aw ( OO; 6;) O (7) . 
Le premier terme représente l’écart de température constant di au signal - 

et le deuxiéme I’écart fluctuant dé au bruit. Le montage approprié (en oppo- - 

sition) des thermocouples donne une tension proportionnelle a ces deux tera 
De tout le spectre de bruit que contient le courant d’entrée seule une bande 


Vol. XI, 1960 Redresseurs synchrones thermiques 93 


étroite autour de w est transmise. Sa largeur équivalente’) est: 4 feca ae Rteres 
ou. t est de l’ordre de 0,5 a 20s, On réalise ainsi facilement des largeurs de 
bande de bruit de 1 4 0,025 c/s. 

Par la suite nous étudierons de plus prés quelques cas pratiques: le montage 
avec thermocouples (§ 3), celui avec lampes A filament en tungsténe (§ 4) et 
celui avec thermistances (§ 5). 

Notons qu'il serait possible de doubler la tension de sortie de ces dispositifs 
en utilisant 4 éléments sensibles. L’égalisation des 4 caractéristiques devient 
alors assez compliquée, tandis que pour deux éléments elle est trés simple, 
sinon superflue. 


3. Systeme a thermocouples 


Dans les caractéristiques de thermocouples on donne généralement direc- 
tement la relation entre la tension du couple et le carré du courant: Vj, = 0 I3y. 
Pour un pont a 2 thermocouples (figure 2) ou les tensions de repos s’annullent 
mais les 2 écarts dus aux termes P, s’ajoutent Ja tension de sortie devient 


Z 


J - 0.5k2 10k2 £ 
controle de V, 10 
BV 
4h : 
12 
2,5 k& 
t —: 
TH2 TH2 
Vs 
[Jee 28Q 
cee ‘ 
IIIA 
(TUSTIN) 
ey abs 
Figure 2 


Redresseur synchrone avec thermocouples. 


3) Cette notion n’est utile que si le spectre de bruit peut étre considéré comme blanc a Vintérieur 
de cette bande. 


a 2 we A 


94 JEAN-PIERRE Borer et Pret CORNAZ ZAMP 


1 i 
Ves OL a hag 5b DT app Tat esr Vit dot cos(d@,t— 6). (8) © 


La résistance de sortie est de quelques Ohms et la résistance d’entrée de quelques — 
dizaines d’Ohms. Pour bien utiliser la puissance du signal disponible il faut 
donc utiliser un transformateur d’entrée pour adapter a l’impédance de cathode — 


ou de plaque d’une lampe. 
Le gain de conversion, définie comme le rapport de la tension continue de 


sortie V, a la tension signal efficace a l’entrée V, .4 sera (pour y = 0) : 
path Vine Wha 
Cae: Veep a ay Tr ett ©) 


R, résistance du fil chauffant; R,, résistance du thermocouple. 

On peut définir d’autres grandeurs utiles qui sont: L’impédance de conver-— 
sion, rapport de la tension de sortie au courant signal (I, .,-) et le gain de 
puissance, rapport de la puissance utile de sortie (P,) a la puissance d’entrée (P,) 


V, b - 

og (10) 
Gees Ve Be he f 
BOB AG 2 Rg Te ep eh) OR a eo 


Le tableau 1 donne quelques valeurs numériques (pour un systéme a 4 thermo-_ 
couples les gains sont a multiplier par 2): | 


Tableau 1 + 

SE a a 
r Pure 1 1 } 

ype eo; SF Ta Ge rary on Smears R, v b Gy, Zo Cp 

Philips Q | mA Q approx.s | V/A? - VIA 4 

. 

TH1ouTH91| 75 5 5,5 0,5 120 |0,8 -10-2] 0,6 {1,1 4 
TH 20u TH 92| 23 10 3 0,8 30 |] 1,3: 20-2") 03m |aleeg 1 
TH 3 ou TH 93) (7,3 20 3 15 7,5| 2,05 10-2] 0,15 |1,28 -1 
- 


F ie gains sont proportionnels au courant de référence; on choisira celui-ci 
e fagon a rester dans la partie linéaire de la caractéristique V. = f(I2 
tableau 1). é Gari! aa 
La figure 2 donne un exemple d’un tel redresseur synchrone tel qu'il a été : 
utilisé dans nos recherches en résonance électronique et nucléaire depui 
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plusieurs années [4]. Actuellement, l’entrée du signal par transformateur a été 
remplacée par une attaque en cathodyne- pour pouvoir utiliser le dispositif sans 
rotations de phase dans une gamme de frequence trés étendue (1 c/s 4 20000 c/s). 
Le courant de repos de la lampe qui passe alors par les thermocouples ne géne 
pas; il faut toutefois en tenir compte en déterminant le courant de référence 


~ admissible. Nous avons inséré dans le circuit de sortie un réglage de zéro per- 


mettant d’ajouter une faible tension dont on peut ajuster la valeur et le signe, 
afin de compenser une légére dissymétrie des thermocouples. On aurait pu 
également ajouter a cet effet une résistance en paralléle a l’un des thermo- 
couples. 


4. Systeme a lampe a incandescence 


La résistance du filament d’une lampe a incandescence augmente avec la 
puissance qui y est dissipée. On monte 2 lampes et deux résistances (ou 4 lam- 
pes) en pont de WEATHSTONE (figure 3) dans lequel on fait passer, conformé- 


ment a la figure 3: 


G\e WOE BTN VA 


nef 


se 


eX 


- 
“al 


1. Un courant de mesure J,, qui peut étre continu ou alternatif mais de 
fréquence nettement différente de la fréquence du signal; par la suite nous le 
supposons continu, 

2. le courant du signal J, de fréquence ow, 

3. le courant de référence J, de méme fréquence, 


Figure 3 
Redresseur synchrone avec lampes 4 filament en tungsténe; principe. 


Le pont, en équilibre avec J, = 0 se déséquilibre proportionnellement au 
courant signal parce que la lampe L, dissipe la puissance Py + P, et la lampe 


- L, la puissance P, — P,. La température des filaments est donnée par la for- 


mule (7). 


pt Paste 
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Dans le schéma de la figure 3 on pourrait également incorporer la tension © 
continue V,, dans le circuit du signal 4 condition de mettre l'instrument de me- 
sure entre A et C. 

La figure 4 donne la résistance R en fonction du carré du courant efficace [? 
pour une lampe particuliérement bien adaptée a notre usage: Philips 8099 Z. 
C’est une lampe 4 V 4 ampoule vidée spécialement étudiée pour Vapplication 
comme stabilisatrice d’oscillateurs. I] suffit d’une puissance de 0,9 mW pour 
doubler sa résistance. 


0 10 20 30-10-84? 


Figure 4 


Lampe Philips 8099 Z; la résistance R et sa dérivée k en fonction du carré du courant J{; gain en | 
tension G, d’un redresseur synchrone a deux lampes. R: x 10? Q; k: x 10? Q/A?; Gy: x 1. 


Le point de fonctionnement en l’absence de signal sera fixé par les courants — 
de repos: I,, et I,: 
2 ie ‘ Le eff 2 
B= (ey) + (34). 


Autour de ce point on peut écrire 


} 
; 
Re RM obAG)s oe ae j 


di2 (11) 


pour [=I, ; 


Pour un courant signal 7, = J, sin w ¢ nous avons pour la valeur moyenne de 
A(I?) avec y = 0 [voir (3a)] 


ON Le L, —- 1 ey 
be cto rege cer ee (12 
La tension de sortie du pont & 2 lampes devient 


: 

1 

Le k : 
og AT lee J 
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Elle apparait avec une résistance interne de sortie R, = Ro. Les gains en tension 
et en puissance et l’impédance de conversion sont: 


V, k 
PWS OE A Tgp oe OS 
: Vs eff 4 Ro In I, _ 
Fe ak | et 2 
Gy == eile, (—= rat.) it 
“4 = entrée 64 Ri : ( 4) 
V, k 
Z, 7 ea ey Vi Dr opp: 


Pour chaque valeur de Jy, donc de A(J6), le gain sera maximum pour J m= 1, off 
et l’expression du gain devient dans ces conditions, 

a k Tie AR 
v max 4 Ry m 2 Ro 


Ey tage a 


c¢ max a m pe 


=~ 
wo 


G 


On voit sur la figure 4 que dans le cas de la lampe 8099 Z ce gain passe par un 
maximum trés plat pour J,=3,6mA, dot I,, = I, 4 =5,1mA. Le gain 
maximum est de 


G 


gain Ors Le ny CL ZO VIA. (Gy maz = 0,0156 . 


¢ max 


La constante de temps vaut environ 1 s. 


Tableau 2 
ee Ro t lee GC. 2s Gy 
approx. s Q/A2 mA VIA 


8099 Z 85 if 57a LOS srl 0,25 ZU2 5) A560 
222.V¥ 7 W 850 2 BLS = Lop 3,15 OWS eZ) 0/56. 10-*} 


Le maximum du gain est une caractéristique unique et trés précieuse de ce 

type de détecteur: le gain sera en premiére approximation indépendant de 
petites variations de J,, ou I,. (De faibles variations dans I, peuvent en effet 
étre introduites par le déphaseur fixant la phase de référence.) 
La faible impédance de cette lampe (environ 85 Q) rend souhaitable une 
adaptation et il serait trés avantageux d’utiliser des lampes de résistances plus 
élevées. Nous avons fait des essais avec des lampes d’éclairage 115 V, 3 W et 
220 V, 7 W (résistance environ 1000 Q au point de repos optimum). La sensibilité 
en courant du dispositif était accrue d’un facteur de 6 (voir tableau 2), mais ces 
lampes ont des inconvénients qui rendent leur utilisation illusoire: 
1. Les filaments sont supportés d’une fagon non rigide en plusieurs points; 
de faibles vibrations suffisent 4 les déplacer ce qui modifie leur résistance: Les 
fluctuations du zéro du pont atteignent plusieurs millivolts. 
id 
“ZAMP X1/7 


; 


4 


ee 


_ 
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2. Plusieurs de ces lampes sont 4 atmosphére gazeuse. Dans ces conditions, 
le temps nécessaire pour atteindre un équilibre thermique aprés enclenchement — 
est de prés d’une heure et la réponse a un signal d’enveloppe carrée se compose 
d’une exponentielle rapide ayant une constante de temps de 2s et dung 
exponentielle lente de plusieurs minutes (trainage). 

La figure 5 donne un exemple d’un redresseur synchrone avec la lampe 
Philips 8099 Z pour une gamme de fréquence de 50 a 200 ke/s. 


Figure 5 
Redresseur synchrone avec lampes 8099 Z. 


; 


Le courant continu J,, est le courant de cathode d’une triode. La sortie de 
la tension continue se trouve au point milieu du transformateur de la référence, 
ce qui diminue beaucoup la composante alternative. Le reste est filtré par ded 
circuits RC. Ce schéma a été congu initialement pour des lampes a résistance 
de 860 Q (lampes d’éclairage 220 V,7 W) qui ont df étre rejetées pour le 
raisons déja citées. Avec les lampes a 85 Q (Philips 8099 Z) une attaque par 
transformateur serait préférable. Néanmoins le schéma de la figure 5 est 


agréable par sa simplicité. Il est 70 fois plus sensible en courant que le systém 
a thermocouple. 


Lorsque le pont est utilisé avec des lampes amplificatrices, il est important 


de Visoler thermiquement si l’on veut avoir une tré ili 
v e trés grande stabilité 
(dérive de 10 uV/h au maximum) : a 
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Remarque. Au lieu de mesurer la température du filament par sa propre 
résistance, on pourrait utiliser des lampes a 2 éléments, l’un de chauffage, 
_Vautre de mesure. On arrive au méme schéma qu’avec les thermistances A 
| chauffage indirect (figure 6 et §5). Les circuits de mesure et de chauffage 
sont séparés et on peut choisir leurs résistances judicieusement en fonction du 
probléme a résoudre. Nous croyons que des lampes a deux filaments métalliques 
seraient préférables aux thermistances A chauffage indirect. A sensibilité 
comparable ces lampes seraient moins sensibles a des variations de la tempéra- 
ture ambiante et a des variations de J,. De telles tubes existent dans le com- 
merce (chronotron Bendix) avec des constantes de temps entre 6 et 120s. 


Figure 6 
Redresseur synchrone avec thermistances a chauffage indirect ou lampes a deux filaments metal- 
| liques. 


5. Utilisation des thermistances 


Les thermistances (par exemple Philips NTC) ont un coefficient de tempé- 
-rature négatif de — 4°%/°C soit environ 10 fois supérieur en valeur absolue a 
celui du tungsténe. L’élévation de température pour une puissance donné est 
‘pourtant 10 fois plus petite que pour le filament des lampes, de sorte que la 
“sensibilité est du méme ordre de grandeur. L’emploi de ces éléments serait tres 
‘intéressant A cause des valeurs de résistances disponibles et du faible encom- 
-brement. 

Nous avons fait des essais avec les éléments NTC miniatures B 8320.03 P2k2 
ayant une résistance de 2,2kQ a 25°. Ces essais nous ont conduit a rejeter 
_provisoirement leur emploi pour les raisons suivantes: wrvrtd é 

1. Le coefficient de température élevé rend le systéme sensible a de trés 
faibles gradients de température et méme a des variations de température 


j 
e 
c 


A 
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ambiante si les deux caractéristiques des thermistances n’ont pas été égalisées 
avec beaucoup de soin. Il en résulte une instabilité du zéro de quelques diziémes 
de millivolts. 

2. Comme pour les lampes a remplissage de gaz, le temps nécessaire pour 
atteindre un équilibre thermique et avoir un zéro stable, est de Vordre d’une 
heure. La réponse comporte également une partie rapide et une partie lente. Un 
montage approprié sous vide peut certainement améliorer cet état de choses. 
Le type B 832004 P par exemple doit donner des résultats plus satisfaisants. 

Finalement on peut utiliser des résistances NTC a chauffage indirect 
(figure 6). Ces éléments sont montés sous vide ce qui élimine les difficultés 
mentionnées. Ce montage a donné entiére satisfaction. Avec les valeurs sui- 
vantes: 


Elément NTC B 8 320 15 P/330 k: 
Ro =1000, R,=1kO; Vi = 28.V, Very =—1V; 


la constante de temps est de 20s et la sensibilité en courant, Z,, de 10 vig 
avec une instabilité du zéro de + 10 pV. 

La sensibilité peut étre accrue dans de grandes proportions, au détriment 
de la stabilité du zéro, en augmentant V,. 


6. Conclusions 


Nous avons décrit un dispositif détecteur de phase dont la réponse 4 un 
signal sinusoidal J, est du type f(I,) I, cos m, ou J, est un courant de référence 
et p le déphasage entre I, et I,. Il est basé sur l’échauffement d’une résistance 
parcourue par ces deux courants. Selon l’élément sensible utilisé la fonction 
f(I,) prend les formes suivantes: : 


a) Thermocouples: /(I,) = c I,. Ceci est le cas de redresseurs synchrones 
classiques avec réponse proportionnelle au produit de deux courants. Ce Pon 
est aussi utilisable comme Wattmeétre et cos g-métre. 

b) Lampes a incandescence: f(J,) part du zéro a peu prés linéairement et 
passe ensuite par un maximum. Pour de faibles valeurs de J, on a une réponse 
C I,I,, prés du maximum de f(J,) la réponse est du type C I, (cas idéal pews 
un lock-in). 

c) Thermistances 4 chauffage direct ou indirect: Ici aussi {(I,) part au 
zéro a peu prés linéairement mais n’atteint en général pas de maximum. 


Tous les systémes envisagés ont en commun les avantages suivants: . 
Simplicité et dimensions réduites. 
Une grande sécurité de marche méme lorsqu’ils sont utilisés dans _ 
appareils électroniques complexes. 
. 
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Pas de vieillissement. 


Constante de temps intrinséque comprise entre 0,2 et plusieurs minutes 
selon l’élément choisi. 

Parmi tous les schémas envisagés, il nous semble que celui de la figure 6 
réalisé avec deux filaments métalliques est le plus avantageux. Il permet de 
réaliser l’opération J, J, cos p ou C I, cos g, selon le point de repos choisi. De 
plus il a son circuit de mesure séparé du circuit de chauffage, ce qui supprime 
tout filtre HF. 
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Zusammenfassung 


Die vorliegende Notiz beschreibt verschiedene Typen von phasenempfindlichen 
Gleichrichtern, die alle auf dem Prinzip der Erwarmung von Widerstanden beruhen, 
die gleichzeitig vom Signalstrom J, und von einem Referenzstrom I, der Frequenz w 
durchflossen werden. Besitzt der Signalstrom die gleiche Frequenz w wie der Refe- 
renzstrom, so ist die mittlere Heizleistung um I, I,.cos y hoher, als wenn Signal und 
Referenz verschiedene Frequenz haben. Die Temperaturdifferenz wird in einer 
Briickenschaltung von Thermoelementen, Metallwiderstanden oder Thermistoren 
‘gemessen. 
Die besonderen Vorziige dieser Gleichrichter sind: Einfachheit, Robustheit, 
‘kleine Dimensionen, Sicherheit der Funktion, keine Alterung, inharente Zeitkon- 
‘stante wahlbar zwischen etwa 0,2 und 100 s. Bei den Typen mit Metallwiderstanden 
besteht zudem ein Arbeitspunkt, fiir den die Ausgangsspannung unabhangig von 
der Referenzintensitat I, wird. 


(Regu: le 18 aott 1959.) 
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Orthogonal Edge Polynomials in the Variational Solution — 
of Some Boundary Layer Problems in Elasticity 
By GABRIEL Horvay, Schenectady, N.Y., U.S.A.) 


1. Introduction 


One of the simplest boundary layer problems for harmonic functions is 
illustrated in Figure 1a2). On edges IJ, III, IV of the semi-infinite strip 
homogeneous boundary conditions . 


®,, = Oy, Diy 0 (la) 


are prescribed, on edge I the inhomogeneous condition 
D0, y) = fy) (1b) 


is given. [We may assume, without loss of generality, that the arbitrary 
prescribed function f(y) is either even or odd in y.] Assuming, for instance, that 
f is an even function, and writing ® as a product function : 


©, = fal) Bal) (2a) 

fal) = cos, 5 reels yStare (2b) 

the equation : 
V2d = 0 (3) 

simplifies to 
mn? 

ar 4 Bere 0 (2c), 

and has the solution : 
BAN) ees (2d 


which vanishes at infinity and is unity at x = 0. 


The solution of problem (1), (3) is then obtained by determining the Fouriel 
coefficients A, in the expansion : 


Ky) = DY An fal) (4a) 


» General Electric Co., Metallurgy and Ceramics Research Department. 
) We refer to the problem as a boundary layer problem, because nonhomogeneous boundar 
conditions are prescribed only along the short edge x = 0, and significant variations in the solution 
will occur only in the immediate vicinity of this edge. 
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and writing 


D(x, y) “= D>: A_®, ae > A, fn(Y) £,(%) . (4b) 
ve 
$-G =0 
I as 
‘a “ 
Reet 
Missy Sah 
©-f(y),®,-F(y) 9-8-0 
a Figure 1 b 


Problem of the semi-infinite strip: a for the V2® = 0 equation; b for the V4@ = 0 equation, 
where nonhomogeneous boundary conditions are prescribed on edge + = 0, and homogeneous 
, conditions on the other edges. Edge III is at distance a > 1 from the origin. 


When we wish to treat the biharmonic problem 
V4d=0 (5) 


in a similar way we run into serious difficulties. There must now be prescribed 
two conditions on each edge (Figure 1b). Thus we may write, similarly to 
problem (1), (3), 


on IJ and IV: G=G,=0, (6a) 

‘on III: Di Giai, (6b) 

on I: either, . D.=f(y),,. O. = 0; (6c) 
or P=), D. = E(y)- (6d) ° 


{The general condition on J is a superposition of cases (6c, d). Without loss of 
generality we may restrict ourselves to case (6c), and assume that f(y) is an 
even function.] No longer do there exist now product solutions of the problem. 
One may still write (the even and odd) pseudo-product solutions introduced 
‘by Fave [1] and Parxovirtcu [2]* 4), 


(P= g- ™* (cos —weoty siny, 4), = Z,Ast., 
n Vn 3 yy Vn 


2 (7a, b) 

: @, = e-'"* (siny, y — y tany, cosy, y), =3,5,..., 

where 

Z Yn =% +18, (8a) 
3) Numbers in brackets refer to References, page 114. ee 

_ 4) The expressions (7), subject to (8), identically satisfy (5) and (6a, b). The indices n = 0,1 

are reserved for the roots y = 0 of (8b, c). 

s 


E 
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are the nonvanishing roots’) of 
DA Dies (8b) 


Bese al (8c) 


sin2y,-2y,=90, n= 
However, on taking the real (Re) and imaginary (Im) parts of the functions ©,, 
the product nature disappears, and therefore the coefficients in the expansions 


f(y) = 3} [A, Re + B, Im] ®, (0,9) , 
F(y) = Di [C, Re + D, Im]®,,, (0, 9) 


(9) 


7 


cannot be determined by the ordinary techniques of Fourier series. Although 
there are ways of solving the problem (see, e.g., [3-5]) the labor involved is 
very great. ; 
Because of the important role equation (5) plays in stress analysis (see, e.g., 
[6]) and because of the need for simple approximate solutions, a variational 
method was developed in [7] which produces, at a modest sacrifice in accuracy, 
a bonafide product function behavior ; 
D, = fal) 804) (10) 
for the biharmonic eigenfunctions, and thereby salvages also the applicability 
of the Fourier expansion technique. In the present report we shall throw some 
further light on this method which has proven useful in a great number of stress 
problems [7-17]. We shall also comment on some axisymmetric stress problems 
which are similar to but more complicated than (5), (6). 


As is well known (see, e.g., [6, 7]) solving equation (5) is equivalent to finding 


‘| 
2. Biharmonic Eigenvalue Problem of the Semi-Infinite Strip : 
4 


the function ® which extremizes the Dirichlet integral$) : 
co 1 : 
e 2 
2EU =| de [B2, +202, 4+@?,) dy. ay) 
0 —1 


For the finite rectangle the functions with negative ~, are also needed. 


5) In the problem of the very long rectangle (a > 1) only the roots with positive py are use : 
6) E is the ‘elastic modulus’, U the ‘strain energy’. A more general expression is ; 


In a rigorous solution the particular value of y is immaterial 


] 1 t and so it may be chosen as 0, as it 
in (11). An approximate solution will, of course, show depenc : 


dence on», 


2EU = |] (20) — 2.1 —») (Oyq%yy — O4,)] dy de 
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We write 
D(x, y) = ARTY) eel) 5 (12a) 
()=f (yay. (12b) 
Then Fs 
O= 25a =o/ [Dy 4, A feb) 2” 2 | 
(13a) 

+2 >" dD) An AGatp Gre t+ DD, An Ardy fh) bn &] aX. | 

If we choose /,(y) orthonormal, i.e., 
fil) = Ox1> (12c) ° 


then (13a) reduces to 


of [SY ARP e2 ye SS ae 20" (13b) 


If we may assume that the orthogonality of f,(y) insures near-orthogonality of the 
functions f,,(y) and of the functions f(y), then (13b) reduces to the diagonal sum 


oY 


0=2E0U~6 3S 4} if (ef? + 22) ey? + Yi) g2] dx (3c) 


end the expansion coefficients appear outside of the integral sign. The functions 
glx ) may now be determined, once and for all, from the Euler equation 

gY — 20, 8 + Ge) e=0- (14a) 
- 


One finds the two fundamental solutions 


gla) =e" (cosh, x + sin B, x) ; (15a) 
Bre 
h,(x) = e~*k*® Sih, 2. (15b) 
k 
g.(0)=1, (0) =9, (15a’) 
(0) =0, AK(0) =1 (15b’) 


vy, =a, + if, = root of yt — 2G Pyy" + Gf") = 0: (14b) 
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Problem (6c) is now solved by determining the expansion coefficients 


A,= (hl) (16a) 


ty) = D/ Ach) (16b) 


and writing (12a). A similar solution holds for problem (6d). 
The orthonormal set /,(y) is, of course, not unique. One may, for instance, 
choose the polynomials’) 


2 i 

f, = Nz (1 — y?)?, Ni =3?-5-—, 

=N, (1-9) 2 32.527. 2b 

fg = N3(1—y?)?y, Ni= 32-5 «7, oe 
2 ‘ 13 

aN — 97)? (S14, A ye Np oe 


[they satisfy the boundary conditions (6a)] and find the approximate eigen- 
values | 


Dame 2075 ales ae 
3 = 3-656 + 1-538 7, | 

(17b) 
v4 = 5-259 + 2.0627, | 


as contrasted with the roots 


V3 = 3-749 + 1.3847, 
(18) 


Yo = 2-106 + 1-1257, | 
yy = 5-356 + 11-5527, | 
: ) 


of the rigorous equation (8). One may also choose the trigonometric functions 


fy = Nz cos*> y ‘ N=, 
fy = Ng cos*> y sine y, N2=8, 
(19a) 
= NN auc a . gente 
Is 4 COS"> y (1 : cos y), M=7> 


*) For generating function, RopRIGUES’ formula, see [17]. 
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and then one finds the approximate eigenvalues 

Yo = 2-120 + 1-097 7 , 

V3 = 3-751 + 11-3167, 

Yq = 5-354 + 1-5307 , 


And of course there exist infinitely many further alternatives. 

Which of the two sets, (17) or (19), constitutes a better approximation ? 
There are various ways of judging their goodness. On comparing (17b), (19b) 
with the rigorous eigenvalues, the roots shown in (19b) appear to be better; 
but if one judges the /,(v) functions by the nearness to orthogonality of their 
derivatives, then one finds that, for the polynomial set, the scalar products 


<fe fi? <fe f> 
i ee 20 
> GPP? Td GDEP uo 
have the values 
—0-288 = cos107°, 0-141 = cos82°, (21a) 


while for the trigonometric set they have the values 
—0-316 = cos108-5°, —0-164 = cos99-5°. (21b) 


So in this respect the polynomial set appears to be better. 
In Figure 2a, b we plot the derivatives 


fs 2 fs a, eee 2 
=—y(l—y?*), 4 ae a 
ee’ - (22a) 
i- 2 a ese 10 View 
Sens YU v9) (1 2530), 26 Ng b Gaagee © 


of the first two even functions of the set (17), normalized, for convenience to 


Goal y= (22b) 
and compare them with the (similarly normalized) rigorously orthogonal 
qT; , 2 P, = adh 3 v2 
Mz Ms (23) 
4%, Pa 2 2 yt 
f=y(l-y)a—3y") —t=1—-10y¥+11>y¥ 


My es M, 


Rete RAAT RY 


ny =e 
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= (fr Vie 


b 


Figure 2 


| 
| 
| 
| 


Plots of the first and second derivatives of the orthogonal polynomials f,, f,, as given by (17a). 
They are compared with the rigorously orthogonal polynomials T,, T, and P 2, P4, where T, is 
taken as My fg, Ps as Ng fs. 


functions§) ; in Figures 3a, b we plot the derivatives 


2h __sinny, ft _ _cosny, 
2 2 

(24) 
ci. ae 1a | 
aN, = Sinwy + 3sin2ay, PN, =coszy + 6cos2my 


for the set (19), and compare them with the rigorously orthogonal functions 


T,=—sinny, P= >=cosry, 
f (25) 
LS SZ es, P, 606 2b vn 


These comparisons are not refined enough to distinguish the better of the two 
sets (17), (19), but they do indicate that even with a 20° departure fro 
orthogonality one obtains functions which have the proper shape. 


8) P;, are the Legendre polynomials; m,, are normalization factors. 
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oe 


We Gan 


Puneet 


oa 
“3 Figure 3 


lots of the first and second derivatives of the orthogonal trigonometric functions /,, f,, as given 
by (19a). They are compared with the rigorously orthogonal trigonometric functions T,, Ty and 
P,, Py, where Ty is taken as mo fg, Py as Mz f9. 


q 
110 GABRIEL HorVAY ZAMP 


One might also use the closeness of the solutions ®, by (17) or (19), to the’ 
rigorous solutions of given boundary value problems, as the criterion for 
goodness. (Such a comparison has been carried out only between the poly- 
nomial approximation and the rigorous solution; the agreement obtained in 
the few instances investigated is very satisfactory; see [8-10, 13].) 

Clearly, however, the values of the scalar products (20) constitute the proper 
criterion for the goodness of an approximation. If in a problem, the type (13a) 
integral rigorously reduces to the diagonal type (13c), then we have the 
rigorous solution of the problem on hand. If the deviation from orthogonality 
of the derivatives is large, then we have a poor approximation on hand. The 
goodness that may be ultimately achieved by an approximation is determined, 
of course, by the governing differential equation. For the biharmonic eigen- 
value problem the scalar products of the derivatives can never reach zero, 
no matter what set is used, because the rigorous eigenfunctions (7) are non- 
factorable. However, the scalar products should, in some sense, tend to a 
lower bound as one approaches the best product function approximation to (7). 
It is not clear to the author how to go about finding sequences of ever-better 
product function approximations’). \ 

We again wish to point out that the ease and convenience of the variational 
method, as presented here and in [7-17], are entirely due to the fact that the 
Dirichlet integral of the problem may be written, in approximation, as the 
diagonal sum (13c). In consequence, only one Euler equation has to be con= 
sidered for the functions g,(x), rather than an infinity of them. It is also obvious 
that the method is by no means restricted to the biharmonic eigenvalue 
problem, but may be used in other mildly skew problems. Indeed, its use has. 
been found useful in numerous other equations; see, e.g., [18]. 


3. The Cylinder Problem 


New features appear in the cylinder end problem. Here one is to find ‘stress 
functions’ m and ® which obey 


Pp=0, Pb=0 a 


®) We may refer to the biharmonic eigenvalue problem as ‘mildly skew’, because there exist 
sets f, for which the functions f;, and f; are only mildly nonorthogonal. If the set p, =f, g, is not 
deemed accurate enough a first approximation to the true eigenfunctions ® , then fe ay intro: 
duce improved (x) approximations 2 
K 


(ae 
Pn = i onnih eae bn.n+k 
SE ‘ 


with coefficients 4 +p and functions n,n K(*) determined by the conventional Ritz procedur 
see [7]. But the elegance of the method now disappears. 


. hrniis aici eh  oek Ree 


en 


Vol. XI, 1960 Orthogonal Edge Polynomials Jeli: 
so that the associated stresses 


6,=V29—9,+ (1 —2)O,,+7@ 


0, =V? 9 — 9, + (1-2) ®,,, 
(27) 


v 


= V? YP Pz: (4 2 v) @.. . r®.., “10 (4 —2y Sin UA +) V2p , 
t =—9,+7®,,,—2(1—»)@,-—rVO, 


satisfy homogeneous boundary conditions on the surface 7 = 1, prescribed 
conditions on the end z = 0%): 


Git, i) = ty,(%, 1) =O, (28a, b) 
d,(co, y) = T,,(c0, vy) = 0, (28c, d) 
and either 
6(0; *)i==-prescribed;. +,,(0, 7) = 0, (29a, b) 
or 
o.(0,7)=0, 1,,(9, 7) = prescribed . (29c, c) 


In the above, v, the ‘Poisson ratio’ is a material constant (vy ~ 0:3 for most 
metals"). The expressions (27) identically satisfy the ‘equilibrium equations’, 
and (26a, b) constitute the ‘compatibility conditions’. 

The foregoing problem is equivalent to-extremization of the strain energy, 


co 1 
y= see faz | is. 


0 0 


(a, st er “3 o,)? 2 Oo, O06 + (6, + G4) oO, “te 7:,| Y dr | 


and it is shown in [20] and [22] that the method of product function approxi- 
mation 


PnlZ 7) = &nl2) fat), Pnlz, 7) = Gal2) Fal”) » (31) 


where /,(7), &(7) are polynomials, may still be used. But while the boundary 

10) The expressions (27) without the /’? terms, were first given by SApowsky and STERNBERG 
[19]. With the 72 terms included they were first given by the author and J. A. MrraBa [20]. For the 
differential equation approach, the /? terms may, of course, be omitted from (27) at once. In fact, 
there exists an alternate formulation of the problem (see, e.g., [21]), whereby equations (26a, b) 
are replaced by the single biharmonic equation V4 L = 0; and the stresses are expressed as suitable 
derivatives of the ‘Love function’ L. For a variational approach, via stress functions, the Love 
function method cannot be used, and the Sadowsky-Sternberg expressions must first be modified 
by addition of suitable 72 terms. This is discussed at length in {20] and [22]. 

1) The numerical results, given further below, are all based ony = 0:3. 
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conditions (28a, b) may be satisfied by the functions 9, through incorporation 
of a factor (1 — 7?)? in the polynomials /,(7), they can no longer be satisfied in 
this fashion by the functions ®,. For the polynomials F,,(7) the process of ortho- 
gonalization, and satisfaction of (28a, b), must be carried out simultaneously. 
To this end, the following method was adopted in [20]. For the sake of 
simplicity, the p, terms were altogether omitted from (30), and F, was written as 


ie = Ang ai Ane ie to ae n+ 2 le ? (32) 


where the coefficients A,, are as yet undetermined. Then (30) reduces, in the 
usual way, to EULER’sS equation 


WH,GY — GB G"'+ CG=0, (33) 


where .%, “B,, G are rather complicated integrals over F, and its derivatives. — 


The solution of equation (33) again leads to functions of type (15), and the 
eigenvalue equation again has the form 


Di Been G=0. (34) ; 


The coefficients in F(7) are now determined as follows. One coefficient is — 


prescribed (say, it is chosen as 1), the others are fixed by conditions (28a, b) 
and the requirement that «,, the real part of the root of equation (34), be an 
extremum. This may be called a ‘generalized method’ of orthogonalization. It 
leads to the eigenvalues 


Yo = 2-691 + 1-3357, 


v4 = 5:589 + 2-464 7 , (35a) 


as contrasted with the roots 


Vo = 2°1722 + 11-3627 , 


v4 = 6-060 + 1-638 ¢ , (36b) 


of the rigorous eigenvalue equation 


Ro) 4, 20-% 
Ji(y) ey ime 
It is our aim to show in the next section, that the foregoing, somewhat 
questionable, method of orthogonalization reduces for an eigenvalue problem 
related to the second-order equation 2 © = 0, such as the Legendre polyno- 


mials, to the conventional method. However, for the more complicated 
cylinder problem the method described is just one of several possibilities. 


(36) 


: 


4 


j 
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Indeed, a more suitable method of orthogonalization was discovered subsequent 
to writing the paper in [20]. It is shown in [22] that F(r) may be related by a 
differential equation to the polynomial /,(r). Thus, incorporating in f,(r) a 
factor (1 — 7?) which assures satisfaction of boundary conditions (28a, b), and 
orthogonalizing these functions in the conventional manner one obtains a 
definition of orthogonalization also for the F,(r); the polynomials so obtained 
agree closely, but not exactly with those obtained by extremizing «,,. 


4. Orthogonalization of the Legendre Polynomials 
The Dirichlet integral for the 


Pd =0 (37) 
equation in axisymmetric spherical coordinates is 
s ny 
U=2n/ rar | [22 + 73] sind do 
0 0 (38) 
ee hes 
pe. 2 2 Leave -viks YA = 
ee [e+ 5 3] dt (f =cos#). | 
Writing 
& 1 
P=) 1), (=f (a (39) 
. =1 
_{we know that 
Z @,={7 | 2 (40) 
r ; 


are the appropriate solutions of (37), but we shall ignore this information] and 
imposing on /(£) the boundary condition 


fl) =1, (41) 
38) reduces to : 
aq OU = 8 / UG) Re + (CL — 0) 1) 87] dr = 0 (42) 
2 0 
a leads to the Euler equation 
A re" +2re'-Ag=0 (43a) 
= 
> Az Ce. (43b) 
“A ' 

ZAMP X1/8 


ate. 
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In the spirit of the previous section, the orthogonalization of the polynomials” 
| ful) = Ago + Aga C2 + #04 + Aan (44), 


requires determination of the coefficients A,, (A, + 0), so that f,(1) be 1, and 
that A be an extremum. In Figure 4a, b we carry out the graphical determina- 
tion of A, for n = 2,4, using the same type of plots that were used in [20] for. 
determining the coefficients of the F(7) polynomials and the values of «,. In_ 
the present case the plotting, of course, is superfluous; we recognize that Ai in” 
(43b) is the Rayleigh quotient!) for the Legendre polynomials; its oxic 
zation is equivalent to solving the Legendre equation 


(eo ili yk (45) 

We now recognize that, in the cylinder problem, equation (34) constitutes : 

a generalization of the Rayleigh quotient, y? taking the place of A. But the 
connection between the polynomials F,, and the extrema of y, is no longer a 
direct one as it is for the Legendre polyeeas! ; 
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A,-2x3 


-4 =2 0 CS aa Same 6 
1 1 
A=] (1— C*) dP/dt) dt | | P? dé 
P J 
a5 -1 


Piha to bbe 


“2 aba B -4 15/4 0 A,-Ox1 A,=2x3 4 
i 1 
A=[a—<) (apjaty at | [ Pac 
ey =i 


P(t) =1-—b-—c+b0?+cC4 


Figure 4 


Determination of the coefficients in the polynomials Py = 1— b+ 00%, Py=1—b—c+dO +e ge 
from the requirement that they render the Rayleigh quotient A an extremum. 
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Zusammenfassung 7 


Homogene Randbedingungen sind langs der Kanten y = + 1 und ¥ = a (a> 1) 
eines Rechtecks vorgeschrieben, inhomogene langs der Kante * = 0. Zur Lésung 
der Gleichung 7? ® = 0, und besonders zur Bestimmung ihres Verhaltens in der 
Nahe von v = 0, wird von einer Minimisierung des Energieintegrals 


2EU=f [[Gi,+ 20%, + @] dy dx 


ausgegangen. Wenn ®(%, y) durch J’ A, ©, = Y A, f,(v) g,(¥) approximiert wird, 
wo die (von vornherein gewahlten) /,,(y) orthogonal sind und die Randbedingungen 
an y= + 1 erfiillen, so ist die Funktionenschar f/(y) und die Funktionenschar 
f(y) nahezu orthogonal, und das Integral vereinfacht sich zu ; 


7 
A 
. 
2EU~ SAL [[O, +202 +B dy de. 
Nachdem iiber y integriert wird, kann man die Funktionen g,,(¥) als Lésungen einer 
n 


einfachen Euler-Gleichung erhalten; die Koeffizienten A, kénnen mittels der 
Fourier-Entwicklungs-Formel bestimmt werden, und das Problem ist gelést. In 
dem elastischen Endproblem des Zylinders kann man das Integral U durch Ein- 
fihrung geeigneter Produktfunktionen ®, = F,(r) G,,(z) in gleicher Weise verein- 
fachen; jedoch findet man, dass F,(v) sich nicht in der iiblichen Weise orthogonali- 


sieren lasst ; man benutzt zur Orthogonalisierung eine Verallgemeinerung der Idee 
des Rayleighschen Quotienten. 


(Received: July 23, 1959.) 
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First and Second Order Distributions of a Sine Wave 
of Random Phase plus Gaussian Noise 
By Roy Lerpnrx, China Lake, California, U.S.A.1) 


Introduction 


One of the simplest types of random signal encountered in communications 
work is the sine wave of known amplitude and frequency, and uniformly 
distributed random phase. The first order distribution of this stationary process 
is well known. The second order density [1]?) is a divergent series [2]; however 
the cumulative distribution is known and convergent [2]. 

A modicum of realism is injected into the study of detection of such random 
signals by addition of a statistically independent Gaussian random noise. 
Rice [3, 4] has made an extended study of this problem, deriving the character- 
istic function of the second-order distribution and a power series for the first 
order distribution, and numerous other properties of this signal plus noise, such 
as envelope, false-alarm, and zero-crossing probabilities. However, he has not 
given the second-order distribution itself. It may therefore be of interest to 

write down an explicit formula for the second-order distribution in terms of 
‘rapidly converging Bessel series, and a like formula for the first-order distri- 
bution. The main tools contain classical integration formulae and a generous 
application of Bessel function addition formulae. In view of the relatively 


: simple results, a more elegant derivation should be obtainable, but the author 


was unable to find it. 
1. The Rice Formula 
Let y(¢) be stationary Gaussian noise with auto-covariance function 
a(t) = E [y(t) y(t + 7)], and let 
c x(t) = Pcos(wt+ 6)+y(4), (1) 
where the phase angle § has a uniform distribution and is statistically indepen- 
dent of y(#). We then have 


lu, v) = E [eis stereo) 


=, (P ye +v?+2uvcosw t) exp|— me [a(O) (w? + v2) + 2 a(t) u al} 


1) Michelson Laboratory, U.S. Naval Ordnance Test Station. 
2) Numbers in brackets refer to References, page 126. 
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for the characteristic function (Fourier transform) of the second-order density 


p(x, %9), so that | 
: oe —4UxX,—10% 
P(%1, Xs) = (ame ‘ [ple v)e 1 2du dv, | (3) 
ia i 
p(x) = Be [ vu, 0) e—i4* du | . 
i (4) 
: . ju x vi 
= oie | ieee u) exp|— re a(0) a e-iux dy. | ; 
a . 
2. An Auxiliary Result and the First Order Distribution 


Both the density #(x) and the joint density #(x,, x.) depend on evaluation 
of a class of integrals, now introduced. Let the function F, be defined by 


2 72 
F(a, b,%; 2) = J,(az) exp(— 7 — iz), (1) 
where 7 is an integer, a is real, b, z, are complex, and o? = 0. We wish to obtain 
G,(a, b, 6?) = : F (a, b, o”, t) dt, (2) 
and will proceed bycontour integration. Since the integrand is analytic in the 
finite z-plane, the complex integral around the rectangular contour with ver- 
tices (-M, N, N + 1c, —M + 1c) is zero for all real c + 0. We prove that the 
vertical contributions tend to zero as M + oo, N > oo, so that the open 
contour (— oo, oo) can be displaced to (—~co + 1c, oo + ic) without changing 
the result. 

The contribution from the segment [N, N + 7 c| can be written as the real 
integral 


c 


ry. 2 
V, =i} J, (aN +ayi) exp[— - (N+iy*—ib(N +éy)| dy 
0 


ar ; 2 3) 
=ifJ,(aN +ayi)exp[— 5 (89) + yy + 5] (3) 
0 


x exp {i (bb y— b, N—o® N y)] dy, 
where b = b, + 1 dy. 
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From the classical inequality [5] 


atte 


Leaps ee 


exp(|Imw|) , (4) 


we deduce 


2\ In (N22 + In| /2 o2 5 
LY, | </ perm Te exp||a y|— od ee ye b. N)| dy 


=i ee (N2 + ¢2)l"|/2 exp(— ae + dy N) [exo(lay 


0 


(5) 


2 4/2 
thy+=F) dy . 


For fixed a, b,, o”, c, the integral is finite and the factor 


(N? + 2)" exp(— )+o as N-—-oo 


for fixed o? > 0, c, n, bs. Hence Vy +0 as N + ov. Similarly the contribution 
from [— M + 1c, —M] tends to zero as M + oo. Thus 


/ F(a, b, 0%; 1) dt = q F (a, b, 0? 2) dz = ‘i Fla,bo%st+icdt. (6) 
‘We now choose c to simplify the integration. Let c = — 6,/o? and check that 


Bia, b,0%,t+ic) = J, [a(t “2 + 4 (mi b)] 


(7) 
a b b, — 1 b,)? 
cool $B s Sase. | 
Hence replacing ¢ — b,/o0? by t,, 
2 72 
G,, (a, 0, 67) = exp(— at— a) exp(— “5 S ) dt (8) 
4 
Next we call ([5], p. 394) on the addition formula 
b tab 
yj, (at - *24) = - 5 Fula i | o ) (9) 


Ree ee ON EN ES 
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and the generalized Weber formula 


| tale t) exp(— | dt | 
ie =[(2Y"eol- a) tate): moons | : 
0, modd. 


If termwise integration is valid, we have 


G. (a. 6,08, = (22) p(y — ex) eo 


g=—0O 


In particular, for 7 = 0, a P,b=x, 0? = a(0), we find 


be) =a Cae exp(— Ao alah) 


x Su I5q Gee fy (sam): 


It is the occurrence of both P/a(0) and P?/a(0) in the density formula which — 
makes asymptotic expansions difficult. Fortunately, the J, functions are 
extensively tabulated and the convergence of the above series is reasonably — 
rapid. (It is interesting to compare the result above with the power series” 
expansions derived by Rice.) If P = 0, we have the noise only, and 4 

1 He 
ee ie eee ee z(0)) . 
as it should. : 
To justify our operations above, it is sufficient to prove that B < oe, where 


B -S fit t)| | Jn-n (— +8")| exp(- -S*) at. a3) 

By the inequality (4) and the integral [6] | 
. | 
firemen S)a= (yr (EY, al 


ee 
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we infer that 


atlim 


. pS los exp (12 
nay oe 202 ot o? £2 
iB =< ae : Sa he i aia 
SD) | imi In — mi! exp(— “,-) at 
co 


(15) 


b \|n-m| 
20 


la 
i BACON = la is 2 \({m|+1)/2 |m|+1 
Sr See ea 
m= —0O 


a 


Now split the last sum into the three parts corresponding to m = max (n, 0), 
m <= min (n, 0), and min (mn, 0) < m < max (n, 0). Each of the first two is 
easily seen to converge by the ratio test, and the third is a finite sum. Hence 
(14) and its corollary (12) are valid. 


3. The Second-Order Distribution 


We now apply the integral formula (11) of the previous section to the more 
difficult task of evaluating the second- order density given by (2) and (3) of 
section 1 as 

1 ar. 
P(x, %2) = ree. / [hey +v?+2uvcosu t) 
7 —00 —00 (1) 


: x exp{— 5 [a(0) (u2 + 02) + 2a(t) wo] imu —ixyvl dude. 


: Since the matrices 


ie COS@ T 
and 


|| COS@ T 1 


are mutually commuting, we can simultaneously diagonalize the two quadratic 
forms in uw, v which enter (1). Application of the 2/4 rotation 


2 are ; ae Te (2) 
2 Z 
“yields 
Pl, Xa) = caf [Jo (P Vr ( 1 — cosw t) + v7 (1 + cosw t)) 


; 3 
x exp|— = [w2 (a(0) — a(t)) + v2 (a(0) + a(t) )] 9) 


= ty ("5") —t, (77 | du, dv, . 


eke EEN 
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From the addition formula ((5], p. 361) 
sf (2 4+ 2—2% 2% 0080) = oy Ty(21) Im\%) COS (m 8) 


m= —C 


for real z,, 2, we have, ee 0 = 2/2, 22 = u2 (1 — cosmt), 23 = 0} (1 + cosw 7) 


Ts (P ea (1 — cos t) + v7 (1 + cosw t)) 


; i ft O 
a aay phe OT 
=> Cl)" Ja ee /2 Uy sin-5*} i a (P /2 Vy cos- 3) : 
If termwise integration is justified, we insert (5) into (3) and obtain | 
1 OT 44 =%s ~ 
Pty) = aa 2! (-1) ” Com (PY 2 sin ps aa) a(t) a 


Geom (? 2 cos = é Sea a(O) + a(1)) ‘ 
1 P? (a(0) — a(t) cos» 7) 
2x /a*(0) — a*(r) exp(— Say ae) | 


a(0) (af + 48) — 2 a(t) 4% % 
x exp 2 (aX(0) = ae) 
00 P (¥, — *,) sin “ ) 
et Fam—24 a(0) — a(t) (7) 


& (#, + %,) cos S | Pasi “ 
* Tamar \— ay age) / 8 \ 2 (a(oy = al 
z | jee cos* 
** \2 (a0) + a(x) 
on combining terms and using the relation J,(— 7 z) = (— 1)" I,(z). 


Fortunately, the triple series above is reducible to a double series, which i 
turn is reducible to an ordinary series. We note the addition formula ([5], p. 361) 


[, (Vwi + w} — 2 w, w, Cos 0) cos [F tan-1 (ese) | | 


W, — W, cosé 


Me 
oN 
= 
= 
as 


1 


;(@2) cos] 0. 


l 
3 


For 6 = 2/2, we find 


I, (V2 + w3) cos (x tant) = 1S Ty; .4(®1) Ip;(w) (—1)5 . (9 


jim = OD 
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For fixed q, 7, take 


P (#, — *,) sin oe P (# + 4g) cos ~ 
Wo = — . ke ee AO es e 
2 a(0)—a(c) °° 4 a(0)+ a(t)’ 
and check that 
Ty q—27(0o) Cos ((2 9 — 27) py) = (—Ht DS” (—1)” | 
‘ | (10) 
( P (%, — #4.) aie | b (4, + %) sil | 
x IL m—2q a a(0) a a a(t) f te ae a(0) a a(t) ae cae) 
where 
Qi wi + wi, yo=tan® 
1 
Inserting this in the previous triple series yields the double series 
D= Site vee ar(Qo) cos ((2 9 — 27) po) | Ws) [,(W4) , (11) 
q,r=-cO 
where 
e4 sin P= COs- 
W_ = Ws 


2 (a(0) — a(z)) ” 2 (a(0) + a(z)) © 


If the above series converges absolutely, it can be re-ordered. Replacing the 


_ sum over g and ry by a sum over g and = g — 7, we have 


eres. ow) ® EES Wit 


D= pas T,,,(09) Cos (2 % Po) po» (—1)* E(w) Ln (Wa) - (12) 
n= —co q=—-o f 
Invoking the addition formula (9) with 6 = 2, k = —n, we obtain 
I_,, (Ws + ) = oD T,(ws) L,—n(ws) (—1)*. (13) 
q=-0O 


Since J,(z) = I_,,(z), we have the ordinary series 

D= 3%; I, (3 + Ws) Len(@o) COS (2 1 Yo) 
(14) 
i woe En Ve (Ws Fw 4) Ty n(@o) cos (2 n Po) , 


SS 


P where «4 = 1, ¢, = 2 form 21. 
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Combining the preceding formulae, we have explicitly 


P(x4, Xs) = aT ENCE exp( a(0) AS oe My 2) 
x exp(— Eg may) 
« See (= Sea ey) (15) 
conan tan (eh ay | 


( \/ (4 — %2)° sin? (4 + %2)? cos? o 
X loa \P Y ay) aay? a) fate? Jee 


Note that if P = 0, so that only the noise is present, we have 


it a(0) (7? + 48) — 2 a(t) x % 
P(%1, %2) = 2 = (a2(0) — a®(z))12 exp ( 2 (a®(0) — a(z)) 1: (16} 


which is the known second-order joint distribution of gaussian noise (bi- 
variate normal distribution). The resemblance between the series (12) of 
section 2 for (x) and the series (15) for #(x,, x.) is not surprising. 

In the special, but important case where a(t) = 0 for t + O (white noise) 
the formula (15) sinplifies to 


D(%1, X2) = saa exp(— 2 a(0) 


ie (soy yx? + #5 + 2 4% COS@ | ; 


A larger class of Gaussian noise processes for which the distribution (15) 
simplifies is that defined by a(r) = a(0) e~“ cosf t. 

The operations performed above are valid in case ¢,, Cs, Cs < 00, wher 
Cy, Cg, and cg are defined by 


co 


a= 3, | [Vi fy 41) | | Jom(a a) 


x exp(— st ui) exp(— <4) du, duz , 
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es oe ae PH 3% | Tym =39(@) | | I, means te) | | Z,(c) | | I.(d) | (19) 


M=—CO g=—0O r=—00 


= D'S \teq-2r(a)| [L()| [20 | (20) 
q=-0CO r=—-0CO 
for real ky, kz, 04, Oo, a, b, c, d. 
We have 
co wo | 1 Hy | 12! 2 Uy ||2m| 

ae ts ee — 
NS Bec eee 

Ps exp(— s u,) exp(— has) du, duty (21) 

[A a} lem 
_— ew | 4d ¢2 \leml+ne (2 \lemi+ne 74 (2|m| +1 
7 NE rl Sa: 


from (4) and (14) of section 13. The ratio of successive terms is of the order of 
const (1/m?) + 0, so the series converges. Likewise 


oo os + eyes (Beas 

‘i ee = As i Piet ; 

(5 |q| d\|r| (22) 
~ 2) (2) 


Fel, aed 


a\l2q-27| (b\lal (c\lrl 
esp 2 ta) Aa) a) 


q=-0O r=—0o 


exp(|a| + [5| + le] + la) 


exp(lal+ [alt lel).  @3) 


r. Now split the threedimensional lattice of summation in (23) into the 
18 disjoint subregions determined by choosing signs of m, 2m — 2 q, 2m—2r7, 
4, ry. In each of these, we have a convergent series. es example, when m = 0, 


“4 =m, q =m the essential part is 


zi co co oo am Ba yr 
A Pies Te 2q—2m)!(2r—2m)!r!q! 
_ am parm yrrm 


a =D) de de giant ims Qi trl’ 
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Since 


’ 


(m+q! 5 1 
HG! 


the latter is 
cu T.L (x B y)™ Bt yt 
VID ecg, rere ee | 
(aBy)™ . 
aioe, (m!)2 2 qtrat 2, een Z 


m=0 


the product of three convergent series. 

Similarly, the twodimensional lattice of situated in (23) can be split 
into 6 disjoint subregions determined by the signs of g, 7, 2g — 27, in each of 
which we have a convergent series. This concludes the validation of the whole 
chain of operations by which formula (15) was derived. 

The author wishes to express his gratitude to Dr. D. E. ZirMER of the 
Michelson Laboratory, USNOTS, China Lake, California, for his kind assistance, 
and to Dr. JuLtus SIEKMANN for providing a German summary. 
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Zusammenfassung 


Die Verteilungen erster und zweiter Ordnung einer Sinuswelle fester Amplitud 
und Frequenz sowie mit gleichmassig verteilter Zufallsphase und GauBschem Zu 
fallsgerdusch werden explizit berechnet durch Fourier- -Umkehrung einer Form: 


nach Rice. Wir erhalten als Ergebnis Besselsche Reihen multipliziert mit Gauss- 
Verteilungen. 


(Received: August 1, 1959.) 
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| Heat-Transfer Near the Stagnation Point of a Body 
of Revolution in the Presence of a Magnetic Field 
By Rupotr X. Meyer’), Los Angeles, Calif., U.S.A. 


1. Introduction 


The problem considered here is schematically indicated in Figure 1. An 
electrically conducting, viscous, incompressible fluid flows in the direction of the 
negative z-axis towards a flat plate located at z = 0. The motion is modified 
by the presence of a magnetic field, which can be thought of as being produced 
in part by a suitable current distribution outside the flow field (i.e. at z < 0) 
and in part by the currents induced in the fluid itself. 


Figure 1 


Streamlines and magnetic field lines (schematic). 


It is assumed that the motion is steady and rotationally symmetric with 
the z-axis as axis of symmetry. 
The magnetic field considered in this paper has the property that its field 
strength vanishes as z > + oo. Furthermore, as is known from the correspond- 
ing non-magnetic case, the viscous forces are negligible compared with inertial 
forces at a large distance from the flat plate. Therefore, in this region, the flow 
_is approximated by the well-known stagnation point flow solution of an invis- 
cous fluid. As will appear in the course of the analysis, the origin of the invis- 
cous, non-magnetic stagnation point flow is displaced in the positive direction 
of the z-axis by an amount which depends upon the magnetic field strength. 
This solution serves as a boundary condition to the actual flow at z= oo 
_ The temperature T,, of the fluid at z= oo is assumed to be uniform; 
“similarly, it is assumed that the plate be kept at a uniform temperature 7. 
The properties of the atch ‘i.e. its density 0, electric conductivity o, magnetic 


1) Space Technology Ee iaiatee, Inc. 
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unit mass c are assumed to be constant. 

The principle aim of this paper is to determine as a function of the radius 7 
the rate of heat transfer from the fluid to the plate. The distributions of veloc- 
ity, magnetic field strength and temperature are also found and explicitly 
computed for some values of the characteristic parameters. The rationalized 


mks system of units is employed. 


permeability u, viscosity 7, thermal conductivity x and the specific heat : 
; 
2. Maxwell’s Equations and Momentum Equation 


MAXwELL’s equations for steady state, 


curllH =j, (1) 


curl E = 0, (2) 


divB =0, (3). 


(H magnetic field strenght, B magnetic induction, j current density, E electric 
field strenght) are supplemented by the constitutive equation i 


B=ywH (4) 
(u magnetic permeability) and Ohm’s law for moving isotropic media, 
j=o(E+v~ B) (5) 


(o electric conductivity, v velocity of the fluid). 

Conservation of momentum is expressed by the Navier-Stokes equation for 
an incompressible fluid supplemented by the term jf x B which represents the 
magnetic force per unit volume acting on the fluid (other forces being negligible 
in the considered applications), 


D 
0 Fp = —eradp + V2 v4+jx B (6) 


(@ density, # pressure, 7 viscosity). Together with the continuity equation 
divy = 0 (7) 


these equations suffice to determine, for instance, the velocity and the current 
distribution in the fluid. Since the properties of the fluid are constant, the 
energy equation is not required at this point; it is used later to determine th 
temperature distribution and heat transfer characteristics, taking into account 


convective and conductive heat transfer, Joule heating and viscous dissipation 
in the fluid. 


ay 
C 
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In the present problem, by reason of symmetry, the derivatives with 
respect to the azimuth vanish. Furthermore, 


j=1.= 0, (8a) 

B,= 0, (8b) 

te 0, (8c) 

E,= E7=E,=0, (8d) 


where the subscripts indicate the components in cylindrical coordinates. The 
postulate of symmetry expressed by equations (8) is justified by showing that 
these equations satisfy both the basic equations and the boundary conditions 
of the problem. Equation (8a) means that the current flows in circular loops 
about the z-axis. Equation (8b) follows from it by an application of Stokes 
theorem to equations (1) and (4). It also follows from equation (8a) that the 
magnetic force jx B has no g-component and that therefore the fluid, which 
initially (i.e. at z= co) has no g-component of velocity, will never acquire 
one. Furthermore, since in Ohm’s law, equation (5), 7 and vx B have only 
gy-components, it follows that consistent with equation (8a) E can have at 


_ most a g-component. But even E£, must vanish, as follows from equation (2), 
which holds at all points. 


Eliminating the current density from equations (1) and (6) by means of 
Ohm’s law, 
curlB = uo (v x B), (9) 


0 oy = —gradp+ 7 V2 y+ . (curl B) x B. (10) 


Solving these two equations together with equations (3) and (7), is equivalent 


to solving the system of equations (1) to (7), since E = 0. 


Written in component form, 


eB BR Sno, 87, 8): (11a) 
Ou, Ou, Op 0*u, Th Oe Uy 070, 
0 (2, or mts a) = +0 Fe % y Or y * a) (11b) 


B, (oe * 2B.) 


i be 0z or 


0?v i” 00; 0*v, 

= Sieh ( : t : a Oz? ) 

0z Ov Y Or Be (116) 
13 jatar (OB, ) 
be ( Oz OV am 


7s 
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the corresponding equations for the other components being satisfied trivially, 
by virtue of the symmetry relations equations (8). Furthermore, 


Bai (11d) 


OU», OP z Ov, As (11e) 


A. Similarity Solution 


An exact solution of the set of equations (11) can be found by substituting 
the trial solution 
B,=7 g'(z) , LepeS —2 (2) ? (12a) 


v,=rf'(2), 1. =—2f@), (12b) 
by — p = hie) + P Ale), (12c) 


where f, g, 2, k are as yet unknown functions of z alone and where the pressure 
is normalized by means of the pressure f, at the origin. Equations (12) are the 
rotationally symmetric analogue to NEURINGER’s and McIrroy’s [1]?) solution 
for the two-dimensional case, and are similar to the non-magnetic case (e.g. 
SCHLICHTING [3]), except that k(z) can be replaced by a constant in the latter 
case. 

Equations (11d) and (11le) are already satisfied by the choice of B and v. 
Substitution into equation (11a) results in an equation in which all terms are 
proportional to 7, yielding therefore, 


gi + 2wo (fg — fg) =0. (13a) 
Similarly, equation (11b) requires that 


gmbani gh 2k 

: ee, 1 
rani 5m (13b) 
where v = n/o is the kinematic viscosity. Equation (11c) yields terms which are 
independent of v and others which are proportional to 72. The coefficients of 
both terms must vanish separately, giving 


y i ay, } fe (f’)2 


20 i +4t i ——h=0 
or integrated 


h—h=2 (+f) (13c) 
and 


ge" —pk'=0 


*) Numbers in brackets refer to References, page 145. 
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or integrated, 


els rs 
pis Malar or | ea (13d) 


where h and & are constants of integration. The trial solution therefore results 
in four equations which the functions /, g, 4, k must satisfy. The manner in 
which the external current loops at z <0 must be distributed, in order to 
produce the magnetic field required by the similarity solution, is not examined 
in this paper. It is noted that it follows from equations (12) that 


jo = 9 (v, B, — v, B, = 207 (f' g—f 2’) 


and therefore the current density in the fluid is proportional to the radius. The 
integrated magnetic field strength due to these currents alone diverges for a 
stagnation-point flow extending to infinity in the radial direction. However, in 
any actual application to bodies of revolution, the radial extension of the 
stagnation-point region is finite, viz. of the order of the radius of curvature of 
the nose of the body. The question of how to distribute the external current 
loops (electromagnetic coils) is therefore coupled with the problem of the flow 
in the large which must include the region outside the vicinity of the stagna- 
tion-point. 


B. Boundary Conditions 


The simplest boundary condition which can be imposed on the magnetic 
field at z = O is to require that B,(0) = By, where By is a given constant. This 
results in a boundary condition on the function g, 


g(0) = — 5 Bo. (14a) 


The radial component B, cannot be independently prescribed; it is obtained as 
a part of the solution. 

Since E = 0 and since there is no current component normal to the plate, 
it is noted that it is immaterial for the solution whether the surface of the body 
be electrically conducting or not. 

Both components of the velocity must vanish on the surface. Consequently, 


f(0) =0, (14b) 
(0) =0. | (14) 


Furthermore, from the definition of £, as the pressure at the origin, h(0)=0, 
and therefore from equation (13c) _ 
b= OF 
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At z= o the following boundary conditions are imposed: B, is required 


to vanish, yielding the condition 
goo) = 0. (14d) 


For a suitably well-behaved function g(z), the latter condition implies that 
g'(co) = 0 and that therefore also B,(r, z > 00) = 0. 

At a large distance from the plate, the condition is imposed that the flow 
be identical with the inviscous, non-magnetic stagnation-point flow. Therciore 


v,(co) = ar, hence j'(oo) =a, (14e) 


where the constant ‘a’ is determined by the characteristics of the flow outside 
the stagnation point region. It follows by integration (cf. the Appendix for the 
asymptotic behavior of the functions / and g), 


v,(co) = — 2a(z—2), 
which is the correct expression for the inviscous, non-magnetic stagnation- 
point flow, and where the constant of integration z = lim [z — /(z)/a] is the 
2—>0CO 


amount by which this flow is displaced in the positive direction of the z-axis 
(due to the action of viscosity and magnetic field). ° 


It follows from equation (13d) that (oo) =k. The product f(z) f"(z) in 
equation (13b) tends to zero as z > oo and therefore 


an ee 
k= 7 4 


from this equation. Consequently, the expression for the pressure assumes for 
large values of z the form 
b+ 5 (GF +07) =p) —2na. 


The flow is seen to satisfy asymptotically the Bernoulli equation, with a total 
pressure which is smaller than the pressure at the origin by an amount 2 n a, 


just as in the non-magnetic case. 
C. Summary of Equations 


Eliminating &(z) from equations (13b) and (13d), and introducing the dimen- 
sionless variables 


C=~ucaz, (15a) 
F(¢) = /=* fe), (15b) 
2 


1S) Sie 2) (15¢c) 
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one obtains two equations for F(¢) and G(o), 
Sty SAT oan (16a) 


m 1 , > BH , 1 " 0 
ge Pe AEP el ee 1G GS (C= Oy (1b) 


where the prime now indicates differentiation with respect to ¢ and where the 
dimensionless constants « and # aie defined by 


a 2 Rn 
Ser 1) OMT Ee re 
hegitice. (17b) 


The quantity « represents the ratio of the viscous shear stress 7 a to the 
magnetic pressure Bj/2 yu. It is also expressed in terms of the magnetic 
Reynolds number R,, = a R? uo (where R is an arbitrary reference length) 
and the Hartmann number H = B, R Vo/n. 

As usual, the quantity f is introduced, which depends only on the proper- 
ties of the fluid. It represents the ratio of the viscous diffusivity » = y/o to the 
magnetic diffusivity (u o)7}. 

The boundary conditions in terms of the new variables are 


alee Ole ie ae Fe () 0 Gt 


; (18) 
BLSE-= cont Sepily eG. 0. 
3. Equation of Energy 
The equation of energy is of the form 
he te 
a Nl Meters (19) 


(c specific heat per unit mass, T temperature, x thermal conductivity). The 
terms on the left represent the rate of heat removed per unit volume by 
convection and by heat conduction, respectively. The terms on the right 
represent the rate of heat added per unit volume through Joule heating and 
through viscous dissipation, where the dissipation function @ in terms of the 
shear stress components 7;; is given by 


dv, - Ov, Ov; 
Dirty eed ty = (ae + oe) 
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In cylindrical coordinates, making use of the symmetry relations equations (8), 


OL oT 0?T a Owe 02?T 
oc (2, ar +4, 7) x ( oy? Se Fimoer te ai on) 
\ 


—Ligenfe[(Se) + C+ GAT + e+ 


This equation is satisfied by an expression which is analogous to eaviation 
(12c) for the pressure, 


(20) 


Ty — T = U(2) + 7? mz) , (21) 


where Jy is the temperature at the origin of the coordinate system. 

Noting that from equations (1) and (12a), 7, = (u)~1 7 g" (2), it is found that 
substitution of the expression for the temperature into the energy equation 
results in an equation which contains terms which are independent of 7, and 
others which are proportional to 7, resulting therefore in two equations for 
1(z) and m/(z). Introducing the dimensionless quantities 


L(¢) = ete (2) , (22a) 
M(¢) =+—> m(z), z (22b) 


where 7, is the temperature at a large distance from the plate, one obtains 
after some calculations, 


1 N" 2 , , fr’ 1 N 
Da, 26 / 4 F! 


where Pr is the Prandtl number, Pr = cy/x, and where the dimensionless 
constant 


d=c(T,,—T) =". 


The first term on the right of equation (23a) arises from the viscous dissi- 
pation, the second from Joule heating. Since these terms are known from the 
previous calculation, equation (23a) is a linear equation for the determination 
of M(¢). Similarly, the last equation then permits to determine ENG). 

If it is assumed that the plate be kept at a constant temperature T = Ty, 


it follows that the fluid also assumes this temperature, as z > 0, and that 


therefore, ; 
at C=0: L=M=_ 0. (24a) 
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The temperature, T,., of the fluid streaming towards the plate is assumed 
to be uniform at a large distance from the plate. Consequently, 


at f=co: L= +i M=0. (24b) 


4. Boundary Layer Approximation 


In almost all foreseeable engineering applications of magnetohydrodynamics 
to fluid mechanics or aerodynamics, 8 < 1. In the case of liquid metals, f is of 
the order of 10~*. In a typical application to the ionized air behind a normal 
shock at hypersonic velocity, B is also of this order of magnitude. Only in the 
case of extremely low densities, or very high temperatures approaching thermo- 
nuclear ignition temperatures, can f be of order one or larger. 

It is convenient to introduce the ratio 


ee eee teak 25) 


which can also be written as y = H?/4 R,. R,, is the viscous Reynolds number, 
R, = a R?/v, where R is a characteristic dimension of the blunt body. Since H? 
represents the ratio of magnetic to viscous forces, and R, the ratio of inertial 
to viscous forces, y represents the ratio of the magnetic forces to the inertial 
forces. ji 

In general, the magnitude of y is limited by the magnetic field strength 
which must be provided and which, in practice is subject to limitations. On 
the other hand, in order to significantly alter the flow as it would be estab- 
lished in the case of no magnetic field, it is in general necessary that the 
magnetic forces be at least comparable in magnitude to the inertial forces. The 
magnetohydrodynamically important case is therefore the one in which y is of 
order one. 


A. Velocity and Magnetic Field Distributions 
Since « and f are both small, one is led to a singular perturbation problem. 


Outside a thin boundary layer, in which the functions undergo rapid changes, 
the first term in equation (16b) can be neglected with the result that 


Gi +2(R Gl-F G)=0, (26a) 
2F BE" — (Fi)? +1—y (26, G — (G)) =90, (26b) 


where F, and G, are functions of ¢,, and where the subscript refers to the region 
outside the boundary layer. In this region, F = Kh, @=G, and €=C,. The 
neglected term represents the effect of the viscous stress, which indicates that 
outside the boundary layer the only forces present are inertial and magnetic 
forces. 
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Since the boundary condition F’(0) = 0, which expresses the vanishing at 
the plate of the radial velocity, is no longer applicable if the viscous stresses 
are absent, one has the boundary conditions 

Atecae= Oc 1 2au0h, G,=+1; 
1 1 1 (27) 
¢; =o: E=-+1, Gi = Ue 


The validity of the boundary conditions at ¢, = 0 is established in che next 
section. The solution depends upon a single parameter, y, which is of order one. 
Functional values such as Fi’ (€, = 0) depend upon y, through the solution, and 
can be expected to be of order one, a result which is of importance for the 
discussion of orders of magnitudes, which follows. 

If the quantity 6 indicates the order of magnitude of the bone layer 
thickness, made dimensionless in the same manner as ¢, the following estimates 
are obtained for the functions F and G and their derivatives, inside the boundary 
layer: 


(ye == 00h, (e) G=41, 
(b) F’ = 0(), HOS EO. 
() F’=0(5), ie "GP BOW 


(4) F" = Ole) 


The justification for the orders of magnitudes (a) to (d) is the same as in the 
case of non-magnetic boundary layers and, therefore, does not require any 
further discussion. 

G is + 1 on the plate and of order one at the edge of the boundary layer. 
(Actually G = + 1, as shown by the considerations which follow.) Therefore, G’ 
would be expected to be at most of order (6). Substitution of these orders of 
magnitudes into equation (16a) then yields equation (g). From this it follows 
by integration that inside the boundary layer, G’ can be approximated by a 
constant, the magnitude of which is determined by the inviscous solution 
[equation (f)]. A further integration then leads to the conclusion that G can 


be approximated by + 1 inside the boundary layer [equation (e (e)]. Finally, 
substitution into equation (16b) gives 


= lx) = o(yB) 
Equation (16a) is seen to reduce to 
G"— 2 F’=0 


inside the boundary layer. With this equation and with equations (e) and (f) 
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above, the function G can be eliminated from equation (16b). With the intro- 
duction of 


sp= po", (28a) 
h=pur, (28b) 
where the subscript 2 refers to the interior of the boundary layer, the equation 
et Gr, Bi— ) - 4y F = —p (cy 1 (29) 


The right-hand side is known from the solution of the system of equations (26). 

The origin of the various terms in this equation indicates that inside the bound- 

ary layer all three types of forces — inertial, viscous and magnetic — play a role. 
The boundary conditions are 


at fos 0: Fe, SO. (30a) 


At ¢,= oe one has the condition Fy(¢,) = F/(0). It can be replaced by the 
weaker condition at {j= 00: BR” =0, (30b) 


because it follows from this condition and from equation (29), 
(Ey +4yF=7(GO)P+1 for ¢,-00, 


since Ff, F’ tends to zero with increasing argument. On the other hand, it 
follows from equations (26) together with the boundary conditions at ¢, = 0 
that 

AGi(0))? +1 = (0) +47 RO 
and therefore 


(EP +4y KR =(HOP+4y RO) for 0. 
The solution of this quadratic equation which applies is 
(6) =H(0) for 2-00, 
as was to be shown. 
B. Temperature Distribution 


Since |x| <1 and |f| <1, equation (23a) outside the boundary layer can 


_ be simplified to 
F, My — FM, = 3(Gi) ee 


in a notation which is self-explanatory. Since the order of this equation has 
been reduced to one, only one boundary condition can be satisfied, viz. 


at f,=00: M,=0. (32) 
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As a consequence of the singular perturbation, M (¢) drops from a value of 
order one at the outer edge of a boundary layer to zero at the wall. It is well 
known from the non-magnetic case, that the thickness of this thermal boundary — 
layer is of the same order as the one of the velocity boundary layer, assuming ~ 
that the Prandtl number of the fluid is of order one. 

This equation can be integrated, in terms of a quadrature, yielding 


co 
« OB Me ts 
M,= Gem (=) a, Te, 
where C is the constant of integration. It follows from the asymptotic behavior — 
of the functions F,(¢,) and G,(¢,) that the integral exists and that the boundary 
condition equation (32) is satisfied, provided C = 0. 

In particular, for €¢, = 0, one has 


M,(0) = —7 Jim | F(6,) i (RY as ; 
a 


As the lower limit of integration tends to zero, the integral becomes unbounded. 
However, the above limit exists, as shown by an application of the Bernoulli- 
L’ Hospital rule, from which one obtains 


M,0) = - 3. 


It follows from equation (26a) that this result can also be written in the final 
form 
M,(0) = —2y Fi(0) . (33) 


By an analogous consideration, it follows from equation (23b) that 
L,(¢,) = 0. It then follows from the boundary condition at infinity, 


independent of ¢,. 


I nside the boundary layer, the following estimates are obtained for the 
magnitudes of the functions L and M and their derivatives: 


(a) L =O(1), (d) M =O(), 
(by = o(-;) (e:) M’ = o(5) ; 
ait (5) (f) M" = O( se) 


The justification for these estimates is analogous to the one given previously 


for the estimates of the functions F and G, and does not differ from the non- 
magnetic case. 


Vol. XI, 1960 Heat-Transfer Near the Stagnation Point 139 


Substituting into equation (23a), it is found that all terms in this equation 
have the same order of magnitude and must be retained. Therefore, Joule 
heating, viscous dissipation, convective and conductive heat transfer are all 
of importance in the magnetohydrodynamic boundary layer. Putting 


M,=M, (35) 


this equation becomes, 
1 " b , 2 , 
5, Mi +2(R Mi— F M,) = (BY = 4 (R) (36) 


if use is made of the equation preceding equation (28a). The boundary con- 
ditions are 


atte Oy Mes () (37a) 


and M,(¢,) = M,(0) as ¢, > co. However, this condition can be relaxed to the 
requirement that 


at t,=00: Mi=0 (37b) 


because if this latter condition is substituted into equation (36) [together with 
results obtained from the examination of the asymptotic behavior, viz. 
I, M3 +0, MZ +0, F > F/(0) and A! > 0 as €, > ov] one obtains 


—M,(t.) =2yF(0) for €,—-co 
and therefore, from equation (33) 


M,(C.) = M,(0) for ¢,—-0o, 
as was to be shown. 
Considering next equation (23b) which determines the function JL, it is 
seen from the estimates of the magnitudes of L and its derivatives in the bound- 
ary layer, that both terms on the right can be neglected. Putting 


L,=L, (38) 
one obtains 
Ea 2Prt, l=, (39) 


where the boundary conditions remain unchanged (replacing ¢ by ¢, and 
L by L,). The solution of this equation is, just as in the non-magnetic case, 


fo t 
Doel f dt ae 2 Pr ib F,(s) dst, (40) 
0 0 


where already use has been made of the boundary condition at the origin, 
and where K is a constant of integration. Putting L,(¢,) = K J(¢,) where I(¢,) 
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is the above double integral, it follows from the boundary condition at infinity 


that K = [I(oo)]-! since the integral converges. Furthermore, the following 


result, needed in the determination of the heat transfer to the plate, is obtained, 


from differentiation of equation (40). Therefore 


1 


fa exp - 2 Pr | Fy) a} Gs 


5. Heat-Transfer 


The rate of heat g transferred per unit area from the fluid to the plate is a 
function of the radius and can be written in the usual notation 


g(r) = hr) + (Ty, — To) » (42) 


where / is the heat transfer coefficient. As is common in the literature on heat 
transfer problems, the dimensionless Nusselt number Nw, defined here as 


Nu(r)= h(r) — bes, (43) 


is introduced. Since 


oT 


q(r) = (S-)_, , 


one obtains from equation (21), after a short calculation, 


Nu(r) = Li(0) — Mi(o) —-“ 


6 (Log — Ty) 
Defining the quantities Nw and Nu) by 
Nu = +15(0) , (44a) 
Nu = —M‘(0), (44b) 
one has 
Nu(r) = Ny Gon a ae 
u(r) Ww) +. Ny o (fot): (45) 


It follows from the derivation that Nw and Ny) depend only upon y 
and Pr. The numerical computations are carried out fot Pr = 0-700. 


t 
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6. Numerical Results 


Computations were carried out on a 1103 Remington Rand digital computer. 
In order to determine the heat transfer data, a solution of equations (26) had 
to be determined first. The result, representing the distributions of velocity 
and magnetic field strength outside the boundary layer is given in Figures 2 
and 3. 


ae 


pha TA Ne nT 
Oe ae 
O2 Oe PECOG TO ITA Th FE Fe 20 22. 24 


eee 


Figure 2 
Function F,(f,). 
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Figure 3 
Function G,(&). 
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Next equation (29) was numerically integrated. The result is given in 
Figure 4. The temperature distribution in the boundary layer can then be 


16 = = 


| 
$ 
. 
| 
| | i 
5 
06 
0-4 T 106 
pls A f ae 
us i LAA ae 
0 204 06 08 90 12 FROTE™ TB eee 
Saba 
Figure 4 : 
Function F,(¢). 
determined from equations (36) and (39) yielding also the Nusselt numbers 


Nu and Nu®) (see table). 


Computed Nusselt Numbers Nu and Nu for Pr = 0-70 


| 


y | Nu Nu 


0 0-666 0-394 
0-274 0-662 0-789 
0-476 0-660 1-057 
0-979 0-653 1-650 
1-48 0-645 2-158 
6:39 0-582 5-274 

106 0-350 16-200 

356 0-264 23-900 


It is noted that the heat transfer at the stagnation point itself decreases 
monotonically as the magnetic field strength is increased. This result is in 
qualitative agreement with NEURINGER’s and McILrRoy’s conclusions [1] for 
the two-dimensional case. On the other hand, Nw) increases, indicating an 
increased slope of the heat transfer rate considered as a function of the radius. 
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The non-magnetic case, ie. By = 0 and/or o =0, corresponds to y = 0. 
The numerical values for the velocity distribution (function F,) in this case are 
in excellent agreement with FROssLING’s work [3]. Nw for y = 0 differs from 


results obtained by REsHoTKO and CoHEN [2] for the non-magnetic case by 
less than 0:5%. 


The author is indebted to Dr. Lez O. HEFLINGER of Space Technology Lab- 
oratories, Inc., under whose supervision the digital computer was programmed, 
and who also suggested the asymptotic expansions derived in the Appendix. 
Thanks are due also to Mr. W1it1AM B. Busu for a number of helpful discus- 
sions. 


APPENDIX 


In this section, the asymptotic expansions for the functions F(¢) and G(¢) 
are found. 


F(¢) and G(¢) can be represented by 
F(o) = oF () + 1 (6) + oF (0) +e: (A-1) 


G(?) = 1G(2) + 2G(2) + °°, (A-2) 


where for sufficiently large ¢ 
[oF (2)4 > |sF(C)| > |2F(O)|>--- 


and similarly for the functions ,G(¢). The zero-order approximation to F(¢) 
is taken as 


oF) =C = 6, (A-3) 


where ¢ is a constant. Introducing the expressions (A-1) and (A-2) into 
equations (16a) and (16b), and neglecting small terms of higher order, one 
obtains 


1G" +2(f- 0) ,6'-2,6=0, oN 
fern tig Te 3 " , I a ANE = 7 
FN + FC — 2) aF SS yertees [2a 1G GG. OF (A-5) 


where the argument of all functions is ¢. These equations together with suitable 
boundary conditions at infinity are taken as the (exact) equations defining the 
functions ,F(¢) and ,G(¢). 

Equation (A-4) which determines ,G(¢), as well as the homogenous part of 
equation (A-5) — considering ,F’(¢) as the unknown — are of the form 


Vea Phe =o)’ —Vieo > 0); (A-6) 
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where k is a constant. The general solution of this equation is 
¥() = ye" = Pax [1 -O@)] + 0%, (A-7) 


where « =Vk (ee a) and where y, and y, are constants of integration, P(x) 


is the error function, 
x 


(x) chats dt. 


For large values of x one has the asymptotic expansion 


re 2 1 3 
yx x [1 —D(x)] ~e* E Kur es a Get | 
from which 


2 z 9 
YQ~gnet rt [1p rtt--| tye. (A-8) 


Therefore, replacing y, by the constant ,g. y,. by zero fas required by the 
boundary condition on G(¢) at infinity], and k by + 1, 


1G) ~Z wer’ C-H4[1-FC-H++--]. (AD) 


Similarly, it follows from differentiation of equation (A-7), that the expression 
which occurs in the inhomogeneous part of equation (A-5), can be written 


24G1G" — (6)? ~ gre 28 -F (¢ C21 —2 (6 —H2+---). (0) 


The solution ,/,(¢) of the homogeneous equation corresponding to equation 
(A-5) is obtained by putting k = 1/6 and by noting that y, = 0 in view of the 
boundary condition on F’(¢) at infinity. Consequently 


Fi(C) wate © (¢ — 2-8 [1-58 C—Ge4..], (atl) 


where ,f is a constant. 


A particular solution of the inhomogeneous equation is obtained by con- 
sidering the expansion 


ENO) We 0-H (¢— Bm Sa, (CB, 


4=0 


where n and the coefficients a; (where a) + 0) are determined by substituting 
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the expansion into equation (A-5), The result of this substitution is 


8 ay (2- =) (¢ — cyte 8 a, (2 - 7) te ~ f)ntt 
+ freee ne 5H] rend] ede 


ee eS Se AE ae Se 


Therefore, the expansion above is seen to satisfy equation (A-5), provided 


nm = —4 and 
( Ls 
o2 18” (30 oa =A 
@=a=++-=0, a=+ 16 ;- &@=—- Tir ete. 
16 (2 - ea eet S- 
2 4 a ( a 
Consequently, 


F'(0) =F (0) +B (0) ~e &-9"8 (6 — 2 [v/ eae Hoey ara t: +] 


and 
a ve ; (A-12) 
4 en 2E-F (¢ C)*) re; (eS ae | 
PEG 
where 


3 
Dy == bg s+ = 0, tg SPs bs SP hs etc. 


a é 5 1 
peg a tee) | C= ——, ogee ay — 


The number of undetermined constants in these equations (C > aie 1g) COTre= 
sponds to the number of boundary conditions on F(¢) and G(¢) at the origin. 
From these results, it follows for instance that 


lim [F() F"@)] =, 
a result which was used in the formulation of the boundary conditions. 
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Zusammenfassung 


Die affine Lésung von NEURINGER und McILRoy [1] fiir die ebene, laminare 
Staupunktstrémung einer elektrisch leitenden Fliissigkeit in einem magnetischen 
Feld wird auf den rotationssymmetrischen Fall erweitert. Die gemass dem Joule- 
schen Gesetz umgesetzte Warme sowie die Warmezufuhr durch Reibung werden 
dabei chne Vernachlassigung beriicksichtigt. 

Die Randbedingungen werden neu formuliert, und Losungen werden erhalten 
in der Gestalt von gewohnlichen Differentialgleichungen, die numerisch integriert 
werden. 

Die Gréssenordnung der auftretenden Parameter wird abgeschatzt in Fallen, 
welche physikalisches Interesse haben. In den meisten Fallen ist es mdéglich, eine 
Grenzschicht-Anndherung einzufithren. Die Resultate hangen dann im wesentlichen 
nur von einem einzigen Parameter ab. Die Geschwindigkeitsverteilung, das ma- 
gnetische Feld und der Warmeiibergang an der Kérperoberflache werden berechnet. 


(Received: June 27, 1959.) 


Kurze Mitteilungen — Brief Reports - Communications bréves 


Eine einfache Verallgemeinerung der Rayleigh-Grenzschicht 


Von Ernst BEcKER, Gottingen, Deutschland 


1. Aufgabenstellung 


Eine fiir alle Zeiten ¢ < 0 ruhende, unendlich ausgedehnte, ebene Wandgrenze 
an ein fiir ¢ < 0 ruhendes, unendlich ausgedehntes, ideales Gas mit der Dichte @p, 
der absoluten Temperatur 7) und der Zahigkeit mo, die im tibrigen der absoluten 
Temperatur proportional sein soll. Prandtl-Zahl Py und spezifische Warme c, des 
Gases werden als konstant vorausgesetzt. Zur Zeit ¢ = 0 werde diese Wand in ihrer 
Ebene mit der Geschwindigkeit wu, = u, t"/t? in Bewegung gesetzt (uw, ist eine 
konstante Bezugsgeschwindigkeit, #4) > 0 eine konstante Bezugszeit, 7 = 0). Gleich- 
zeitig soll die Wandtemperatur auf den konstanten Wert T,, gebracht werden und 
durch die Wand Gas mit der Normalgeschwindigkeit v,, = v, ¢}!2/f'2 abgesaugt 


Op = 0) bzw. ausgeblasen (v,, > 0) werden (uv, ist wieder eine konstante Bezugs- 
geschwindigkeit)!). 


2. Losungsverfahren 


Die fiir t > 0 sich an der Wand ausbildende Geschwindigkeits- und Temperatur- 
grenzschicht hangt ausser von # nur von der Koordinate y in wandsenkrechter 
Richtung ab und wird durch folgende Gleichungen bestimmt Lo SIRye 


@ (uy am Uy) = (u Uy) y ’ (1) 


, 1) Zusatz bei der Korrektur: Inzwischen ist eine Arbeit von Kwanc-Tzu YANG erschienen 
(Unsteady Laminar Compressible Boundary Layer on an Infinite Plate With Suction or Injection 
lie Aero-Space Sci. 26, Nr. 10, S. 653-662 [1959]), in der Grenzschichten der hier peteachteten 
ieee Kontrolle von Rechnungen nach dem von- Karman-Pohlhausen-Verfahren herangezogen 


) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 152. 


~F 
a 


ae ee! RY 
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: ‘at 2 
Q(T +9 T= 5 Ty t+ ,), (2) 


Q + (Q%),= 0. (3) 


Hierbei ist wie tiblich vorausgesetzt, dass die Absauge- bzw. Ausblasegeschwindig- 
keit v,, klein gegen die Wandgeschwindigkeit w,, bleibt, da sonst die Grenzschicht- 
vernachlassigungen nicht gelten, besonders der Druck nicht als konstant betrachtet 
werden darf. Direkt nach Bewegungsbeginn ist diese Voraussetzung wegen v,, ~ ¢-1/2 
nicht erfiillt, aber fiir sehr kleine ¢ sind ohnehin auch fiir v,, = 0 die Grenzschicht- 
vernachlassigungen nicht mehr zulassig [1]. 

Die Randbedingungen werden weiter unten formuliert. Diese Gleichungen lassen 
sich durch folgende Transformationen vereinfachen (vgl. [1, 3]): 


Pat, (4) 


y 
y= | 2 ay, | (5) 
¢ Ow 


Die Dichte 9,, an der Wand ist wegen T,, = const ebenfalls konstant. Dann erhalt 
man 


0 Q 0 
= a 6 
0 a - a 0 a 0 
ee tee O¢ 25 | “= if! 
aor / 04 | Oy ~ Ot Oy Oy (7) 
0 


Beim Ubergang von der zweitletzten zur letzten Zeile wurde dabei Gebrauch von 
der Kontinuitatsgleichung (3) und der Randbedingung v(0) = v,, gemacht. Die 
Gleichungen (1) und (2) gehen nun mit der Beziehung u @ = “, @, (Voraussetzung, 
dass die Zahigkeit der absoluten Temperatur proportional ist!) tber in 


Up + Vy Uyr = Vy Uy yrs (8) 
y v 
peak as ee aie ie a (9) 
: ge ie ee 2 ia og Cte 


mit Yap = My/ Ow: , " 
Fiir das Folgende wird zunachst das der Reibungswarme Rechnung tragende 


letzte Glied in Gleichung (9) vernachlassigt; die folgenden Uberlegungen gelten 
daher nur, solange die Mach-Zahlen in der Grenzschicht klein bleiben. Erst zum 
Schluss wird auf das Reibungsglied in (9) eingegangen. Die Gleichungen (8) und (9) 
lassen sich mit dieser Vernachlassigung durch folgenden Ansatz lésen: 


mit 
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Aus (8) und (9) ergibt sich fiir F(”) und O(n): 


F'+2(n—17) F’—4nF=0, (13) 
Q"” 
eas n—rO=0, 14 
pw t2n—7 (14) 
wobei r ein das Absaugen bzw. Ausblasen kennzeichnender Parameter ist, nadmlich 
Uo Hl? 

is yl (15) 

y > 0 bedeutet Ausblasen, vy < 0 Absaugen. Die Randbedingungen sind 
n=0: F=O=1, (16) 
R= co; =O =O). (17) 


Man sieht nun leicht folgendes ein: Es sei die Losung des Problems ohne Absaugung 
(y = 0) bekannt und mit F(n; r= 0) = f(y), O(n; 7 = 0) = B(m) bezeichnet. Diese 
physikalisch nur fiir 7 > 0 sinnvollen Lésungen lassen sich in das Gebiet 7 < 0 
analytisch beliebig weit fortsetzen, weil die Gleichungen (13) und (14) in diesem 
Gebiet fiir beliebiges y singularitatenfrei sind. Die Lésung fiir beliebiges 7 erhalt 
man dann als 


FQ) = OS, (18) 
e(n) = (19) 


denn diese Lésungen befriedigen die Gleichungen (13) und (14) und die Randbe- 
dingungen (15) und (16). 


3. Rayleigh-Grenzschicht mit Absaugen oder Ausblasen 


Als Beispiel wahlen wir m = 0; in diesem Fall wird die Wand zur Zeit t= 0 


ruckartig mit der Geschwindigkeit u) in Bewegung gesetzt (Rayleigh-Grenzschicht 
[1-3]). Es ist dann 


2 ra 2 
Gh See nds 20 
7 / (20) 
P(e el eee (21) 
Yar 


Der Ausdruck (21) fiir O(n) gilt iibrigens fiir alle » > 0. Die Verlaufe von f(y) (Ge- 
schwindigkeitsprofil) und #(7) (Temperaturprofil) unterscheiden sich nur durch 
einen MaSstabsfaktor, so dass wir uns auf die Diskussion des Geschwindigkeits- 
profils beschranken kénnen. 

f(n) ist in Figur 1 fiir — oo < 9 < oo schematisch dargestellt. Den Verlauf von 
F(n) fir y + 0 erhalt man hieraus nach (18) durch Verschieben des Koordinaten- 
nullpunktes nach 7 = —r und Anderung des Ordinatenmafstabes derart, dass die 
Kurve f(y) die verschobene Ordinatenachse im Punkte F = 1 schneidet Mit 
anderen Worten: Die Grenzschichtprofile fiir alle y < 7), wobei 7, eine beliebige 
Konstante ist, sind simtliche in dem Profil fiir y = 7, enthalten. Da die Kurve /(7) 


rae he 
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und damit auch #7) an der Stelle 7 = 0, f = 1 bzw. # = 1 einen Wendepunkt und 
sonst keinen hat, zeigen alle Geschwiudigkeits- und Temperaturprofile fiir 7 > 0 
(Ausblasen) einen Wendepunkt, dagegen sind diese Profile fiir y < 0 (Absaugen) 
wendepunktfrei. Dies entspricht ganz den Verhiltnissen bei der analogen statio- 
naren Grenzschicht (vgl. [2], Seite 219). Die Geschwindigkeits- bzw. Temperatur- 
profile fiir verschiedene Werte von ¢ sind in Figur 2 aufgetragen. 


Fin) flir r>0 fi) 
| 


Fin) fir r<0 


| 
| 
| 
+ 
| 
\ 
| 
\ 
\ 
| 
\ 
| 
i 
| 
I 
if 


70 70 10 : 
fir r>0 furr=0 firr<0 
Figur 1 


Koordinatensystem und Verlauf von F(7) fiir r > 0 (Ausblasen), 7 = 0 und 7 < 0 (Absaugen). 


Fy), 80) 


Figur 2 


Geschwindigkeits- bzw. Temperaturprofile der Rayleigh-Grenzschicht mit Ausblasen (y > 0) oder 
Absaugen (r < 0) bei vernachlassigter Reibungswarme. 


Fiir die Wandschubspannung erhalt man (Index w bezeichnet Gréssen direkt an 
der Wand) * 


Ty = He Ge (22) 


y=0- 


Nach kurzer Rechnung unter Benutzung von (10) und (12) mit » = 0 erhalt man 
hieraus als Wandreibungskoeffizient c,, 


uf lee las Me PO) (23) 
Be SIZ Ong HG ige § T(=F) 


c 
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Anhand des in Figur 1 wiedergegebenen Verlaufs von f(n) ist klar, dass ¢,, mit wach- 
sendem 7 abnimmt. In Figur 3 ist der Faktor |/’() |/f (—7) aufgetragen. Fiir die 


20 = ia 


| 


A 
f'(r)|  \Absaugen <> Aushlasen 


f¢-P) 
-08 -04 0 o4 08 
fi —. 
Figur 3 


Anderung des Wandreibungskoeffizienten der Rayleigh-Grenzschicht durch Absaugen oder Aus- 
blasen [vgl. Formel (23) ]. 


pro Zeit- und Flacheneinheit durch Warmeleitung auf die Wand wbergehende 
Warmemenge q bzw. die entsprechende Nusselt-Zahl Nu = ¢ up t/{Ay - (Tw — To)} 
(4,, = Warmeleitfahigkeit an der Wand) ergibt sich im Sonderfall Pr = 1 


ae 
Nu = vw. me ; (24) 
w 
wahrend fiir vy = 0, aber beliebiges Pr 
rae 
Nig = “. : ee Pye (25) 
Ww 


gilt (vgl. [3]). 
4. Allgemeinere Falle 


Wenn 7 + 0 ist, dann gelten die in Abschnitt 3 hergeleiteten Beziehungen 
weiterhin unverandert fiir das Temperaturprofil, da m in die Gleichung (14) fiir 0 
gar nicht eingeht. Das Geschwindigkeitsprofil ergibt sich aus (13). Fiir 7m = 1 und 
y = 0 erhalt man zum Beispiel (vgl. [2], Seite 176) - 


[e,<) 


2 2 i 
SSE eaeps (14 20) fer ds—ne "|, (26) 
4 
” 
Hieraus erhalt man nach (18) bzw. der in Figur 1 erlauterten Koordinatenverschie- 
bung auch F(n). Da f(n) nach (26) in — co < 9 < oo keinen Wendepunkt besitzt, 
erhalt man auch hier beim Ausblasen keine Geschwindigkeitsprofile mit Wende- 
punkt. 

Die Beriicksichtigung der seither vernachlassigten Reibungswarme in (9) andert 
an den Uberlegungen fiir das Geschwindigkeitsprofil gar nichts. Dagegen erhalt man 
mit den oben benutzten Randbedingungen fiir die Temperatur (konstante Wand- 
und Aussentemperatur) im allgemeinen keine Temperaturprofile mehr, die nur von 


ae 


Le PATER ENV TA, 
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n abhangen («ahnliche Profile»). Solche Profile ergeben sich nur fiir spezielle Rand- 
bedingungen, zum Beispiel dann, wenn man T,, = T, fordert. Die T emperaturvertei- 
lung kann dann durch folgenden Ansatz beschrieben werden: 


(ae Pac Setlllf pe Se 7 
i 0 | i2n 2 Be (7) , (2 ) 
und fiir © folgt aus (9): 
1 Si +7 ~ i) 
py ? +2(n—7)O—-8nO0=-2F%, (28) 


mit @(0) = O(co) = 0. Der einfache Zusammenhang zwischen dem Temperatur- 
profil fiir y = 0 und demjenigen fiir” + 0 geht hierbei allerdings verloren. Fiir 2 = 0 
(Rayleigh-Grenzschicht) erhalt man im Sonderfall Py = 1 auch fiir beliebige, kon- 
stante Wandtemperatur 7), ein einfaches Ergebnis. Als Temperaturverteilung in 
der Grenzschicht ergibt sich hier: 


B= i+ (= T,) P= Te) SF) (29) 
mit 
up 
T,=To+ 3; (30) 


Man bestatigt dies leicht durch Einsetzen in Gleichung (9). Dies ist ibrigens dieselbe 
Abhdangigkeit des Temperaturprofils vom Geschwindigkeitsprofil, wie man sie bei 
stationdaren Grenzschichten ohne Druckgradient hat (Croccosche Beziehung). Falls 
kein Warmeiibergang durch Warmeleitung zwischen Wand und Gas stattfindet, 
das heisst, falls OT /dy = 0 fiir y = 0, folgt aus (29) fiir die Temperaturverteilung : 


= F(2—F). (31) 


In Figur 4 sind die Temperaturprofile nach (31) fiir verschiedene ry aufgetragen, 
F ist hierbei durch Figur 2 gegeben. 


10 
A r+ 
6 el 
0, S08 
ne +10 
On bbe 
Le 
02 T 
Tp 
af 


Figur 4 
Temperaturprofile der Rayleigh-Grenzschicht beim Thermometerproblem. 
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Zum Schluss sei noch folgende Bemerkung angefiigt: Das in Abschnitt 3 behan- 
delte Beispiel x = 0 kann als Sonderfall einer allgemeineren Klasse von ebenen 
Grenzschichtstrémungen aufgefasst werden, die sich dadurch auszeichnen, dass 
die Aussengeschwindigkeit sowie die Wand- und Aussentemperatur nur von #/t¢ 
abhangen (¥ = Koordinate in wandparalleler Richtung) und die Absaugegeschwin- 
digkeit einem Gesetz von der Form v,, = #1” (x/t) geniigt. In diesem Fall lassen 
sich die instationaren Grenzschichtgleichungen von 3 auf 2 unabhangige Variabeln 
reduzieren. Einen tiberblick iiber solche Grenzschichten fiir den Sonderfall v,, = 0, 
wie sie in Stosswellenrohren und adhnlichen Einrichtungen auftreten, gibt die 
Arbeit [3]. 
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Summary 


A simple generalization of the theory of the compressible boundary layer near 
an infinite flat plate to the case with suction or blowing out is given if at the time 
t = 0 the plate is set into motion in its own plane with velocity u, ~7#”. The normal 
velocity at the wall shall vary with time according to uv, ~7—1/?. In that case one 
gets similar boundary layer profiles for all ¢ > 0, which can be reduced to the profiles 
without suction or blowing (v,, = 0) by a simple parallel displacement and stretching 
of the coordinates. As an example the Rayleigh boundary layer (m = 0, u,, = const) 
is discussed. 


(Eingegangen: 21. August 1959.) 


The Back Effect of a Wall on a Jet 


By Harris C. Levey, Nedlands, Australia?) 


Introduction 


_ There is a great deal of interest currently in bodies which are supported by 
jets and stay relatively close to the ground. Although intuitively one seems to 
regard jets as behaving in a parabolic manner, that is, that no back effects are felt 
at the source of the jet, this cannot really be so for subsonic jets since the equations 
governing the flow are elliptic. 

More specifically the intuitive argument runs as follows. If the source of a jet 
of incompressible and inviscid fluid is at distance h from an obstacle and if d is a 
typical jet width, then if d/h is small the typical length scales for longitudinal and 
lateral variations are in the ratio h/d, not too near the source or the obstacle. From 


1) Mathematics Department, The University of Western Australia. 
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the continuity equation the longitudinal and lateral velocity components are in the 
same ratio. This leads to the parabolic form of the governing equations and the 
result that for impingement problems the jet is assumed uniform apart from the 
influence of the obstacle. 

Of course this is really a singular perturbation problem in the parameter d/h 
and the parabolic solution is not valid near the jet source. However, we feel that a 
more interesting question is how do the gross features of the jet, namely momentum 
and mass flux, depend on the ratio h/d? Accordingly, the simplest possible two- 
dimensional model that can be treated by free-streamline theory is chosen, namely 
an incompressible fluid confined in an infinite chamber whose side is parallel to a 
wall distance h away. The chamber side has a slit of width d through which the 
fluid escapes (Figure 1). 


Figure 1 


The physical plane; with axis of symmetry OA, chamber side DE, wall AB, jet boundary DC. 


The solution is carried through in the usual manner by mapping the flow region 
in the plane of the complex velocity potential, w, and in the plane of the Helm- 
holtz potential, 2, on an appropriate region of the plane of an auxiliary complex 
variable ¢. If J and M are the excess momentum flux and mass flux and J, My 
are their corresponding values when h/d is infinite, we are now able to express the 
ratios J/J,,, M/M,, in terms of complete elliptic integrals. 

On closer examination the interesting result emerges that J/J,, and M/M,, 
remain nearly equal to unity not only for large h/d but down to values of h/d =- 1. 
Of course, no account is taken here of laminar or turbulent mixing, for example, so 
the results although of intrinsic interest can be taken only as indicating a possible 
trend in practice. 


The Solution Details 


The stagnation pressure p, is supposed greater than P, the pressure of the fluid 
outside the chamber, so that the fluid emerges from the slit as a jet, the constant 
speed on the jet boundary being U. The jet divides symmetrically at the wall and 
ultimately flows parallel to the wall in a half-jet of width H. The axes chosen are 
shown in Figure 1, and since there are two stagnation points on the axis of sym- 
metry, at y = co and y = —h, it is assumed that there is one point, X, of maximal 


speed on this axis. 
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If g, 0 are the polar velocity components and 


Z2=H%+0y 
we have 
dw ~i0 
eee i 
dz He () 
and 
Q = logs +18. (2) 


We need consider only the half of the flow field in ¥ = 0 by symmetry. The regions 
corresponding to this in the w and 2 planes are mapped onto the positive quadrant 
of the plane of an auxiliary complex variable ¢, such that the flow boundaries in 
both planes map onto the axes Ret= 0, Im#= 0, the points Cz DEN ORAGEE. 
mapping onto ¢=ioo, 0,1, 1, k-4, oo respectively. Note that the flow region in the 
Q-plane is bounded by a semi-infinite pentagon that has a doubly-described reen- 
trant side O.X A with interior angle 2 z at X. The mapping functions are 
; Cth 1—? 
w=7 UH = log —5 (3) 


and 


1+t 1—kt 


where k < 1. (In passing, the value of ¢ at the point X is found to be k—/?2.) The 
branch of the logarithm is chosen that is real for ¢ < 1 and real, and the singulari- 
ties are avoided by passing above them. 

If 2c is the jet-width at a given value of y, the excess momentum flux (per 
unit breadth normal to the plane of flow) and the mass flux are defined respectively 
by 


c 


J=2/(6- P+ ov) ar, M=-—2ofvdv=20UH, 
0 0 


where v is the y-component of the velocity and @ is the density. From momentum 
conservation j is the excess pressure on the wall, so that 


J =2f@ — P) dx = of (U*— gh ae, 
0 0 


where the integration is along the wall y = —h. When h/d is infinite, it may be 
deduced?) that 
ph zm do U? 
of es m+ 2 ) 
and 
Nie ee 
= Te 
so that 
co 
RE Cee” U2 — @ 
i (Caz) f (Set), mymen 
0 


2) H. Lams, Hydrodynamics, 6th ed. (Dover 1945), 


Es 


GN ina Bad at 


> 
S| 


“Avian Weis SO than aS 


Vol. XI, 1960 Kurze Mitteilungen — Brief Reports Communications bréves 155 


and 


We obtain the relations between h, d, H and k by using 


dz dw 
ee STH 
ae Sa 
so that 
ko 
: I Ye ee 
pete ea PSA en ea SA eye 
fh (7) 


where the roots are real and positive when ¢ < 1 and real, and the path of inte- 
gration detours above the singularity at ¢ = 1. Similarly we find that 


oo 


) - (1 + R) je (2 — 1) 
ko 


‘ee (= + 2 
qe 
(6 + 1)? ¢ — 1) 48 (At 2 1)2? (kt — 1)? de. 
The substitution 
sin @ = is 
Es 
reduces the integral in (8) to the complete elliptic integral of the third kind 
IT(k?, k)*), while the substitution 
(=f sind 
in the integral in (7) reduces it, after some manipulation, to 


(i =) {a + : he (1 f) log + HB + ai k*(1— A) + id (K’— By} 


, 


where K, E are the complete elliptic integrals, with modulus k, of the first and 
second kinds respectively, and a dash denotes the same function of the complemen- 
tary modulus k’[= (1 — k?)/?]. Since JJ may also be expressed in terms of E we have 
finally ‘ 


1 1 7 / 
h a (1—k) + 5 (K’ — E’) 
o- 1) : (9) 
it 
1+ — a) (1—A) 
if. a 2 iy 
and 
w+ 2 
13, 
7 a | D Ie: a 


3) P, Byrp and M. Frrepman, Handbook of Elliptic Integrals for Engineers and Physicists 
(Springer-Verlag, Berlin 1954). 


156 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


where 
eee 
D=14Ee =) logs. (12) 


Now when & is small, 
4 


1 2 
ae E= 51+ Of), 


Ko= log a6 o(z? log +) ‘ B’=1+ O(a? log 


so that the case k > 0 corresponds to h/d > oo, and, moreover h/d varies so rapidly 
with k that when h/d = 1, & is still only 0-01 and then J/J. differs from 1 by 
merely 0-004 and M/M, differs from 1 by 0-014. From about h/d = 1 to h/d = oo 
then, good approximations to (10) and (11) are 


a 


m\-1 1 (2+ a) h 
ee (1+ 5) exp| = it d (13) 
and 
ut —1=—(44 nm) (24+ 2) exp E et See sc =] (14) 


At the other extreme & = 1 gives h/d = 0 and there M/M,, = 0, while 


if ik 
“— = — + 7z1= 0-818, 
Jo 2 


so that there is not a great deal of variation of jet momentum over the whole range. 

The curves are shown in Figure 2. As remarked on before the interesting fact is 
that, as far as gross features are concerned, the jet behaves as if governed by 
parabolic equations down to such comparatively small values of h/d. It seems 
to the author that this result is by no means obvious. 
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Variation of excess-momentum flux and mass flux with wall spacing. 
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Zusammenfassung 


Ein ebener Fliissigkeitsstrahl von-der Weite d wird normal gegen eine Wand 
geleitet, die von der Quelle h Einheiten entfernt ist. Man nimmt gewohnlich an, 
dass, wenn das Verhiltnis h/d gross ist, die Auswirkung des Hindernisses nicht zur 
Quelle des Strahles zuriickwirkt, so dass im besonderen wichtige Bigenschaften des 
Strahls, wie sein Impuls und die Massenstromung, durch den Widerstand unberiihrt 
bleiben. 

An einem einfachen Modell, in dem der Strahl, bestehend aus einer inkompres- 
siblen und reibungslosen Fliissigkeit, von einem Schlitz in der Wand einer unend- 
lichen Kammer ausstrémt, zeigen wir, dass, soweit diese Haupteigenschaften be- 
troffen sind, der Strahl durch das Hindernis beinahe ungestért bleibt, selbst fiir 
Werte der Verhaltniszahl h/d, die nahe an der Einheit liegen. 


(Received: July 29, 1959.) 


Bemerkung zu einem Hauptachsenproblem in der Plastizitatstheorie 


Von Hans ZIEGLER, ETH, Ziirich 


1. Einleitung 


Die Fliessbedingung eines isotropen idealplastischen Materials!), das unter 
hydrostatischem Druck nicht fliesst, ist bekanntlich [2] von der Form 


F(Ss, Ss) = 0, (1.1) 
wobei 
u 2 a 2 at ines 
Se 3 (oz fe aes Oy Oz — -*:) ae si 
1 lied, 
S3 = 27 (2 Oy Oy G,) (2 dy = 0, — G,) (2 C= Og= G5) ( 2) 
1 
3 (2 6,0, > Gyotey —~:+ + 2 Ty 5 Tee Cry 


die Grundvarianten des Spannungsdeviators sind [3] und die Funktion F so gewahlt 
ist, dass F(0, 0) < 0 gilt. So hat im Falle der Fliessbedingung von v. MisEs [4] 
(1.1) die einfache Form 


: Oo. 10K (Ges) 


F=S,-3 


wenn o, die Fliessgrenze im linearen Spannungszustand bezeichnet; nach der 
Fliessbedingung von Tresca [5] ist dagegen 


F = 4 Si — 27 S2 — 9 03 S3 + 6 of S, — of = 0. (1.4) 


Im Rahmen der Theorie des plastischen Potentials [6] lautet das zugehorige 
Fliessgesetz, in den plastischen Verzerrungsgeschwindigkeiten é,..., Yyz) +--+ Se 


1) Beziiglich der Definition siehe [1]. Die Ziffern in eckigen Klammern verweisen auf das 
Literaturverzeichnis, Seite 163. 


158 Kurze Mitteilungen — Brief Reports — Communications bréves ZAMP 


schrieben, wiederum allgemein 


Lines On ee ie 

Pome ae tN OS, Cone ao getoo, os ae 
; 4 OF a(se-- get ee 

Me” Oxy: OSs. , Ot elo aD Semone te 


Die hier auftretenden partiellen Ableitungen von S, und S; sind nach (1.2) durch 


is) ih 
aan =i (2.en—6,—= 0,) umes 
OS, 
= 2 ’ ’ 
OT y« "ua 
IL 
pee eee ge o2 — o2 + 40, 6, — 26, 6, — 2 Gy G,) [ ee, 
00, Oe ee i ; 
lea, (2 a2, = 2a Sy) eee 
0s 2 
ae =->5 (2 Og — Oy —10;) hyp ses De Fcc me | 
gegeben. 
2. Problemstellung 
Ist die z-Achse Hauptachse des Spannungszustandes, dann gilt 
Ty, = Tn = 9, (Zes) 


und aus den Beziehungen (1.6) folgt, dass die partiellen Ableitungen von S, und S, 
nach t,,, Sowie T,,, verschwinden. Somit ist nach (1.5) 


OF OS, OF 0S, ' 


(2.2) 


CHE” GIS OEY OS n\0o 0 
lee “Dee ge OSE. = ; 
Zwar muss grundsatzlich auch die Méglichkeit 2 = co im Auge behalten werden, 
wie etwa beim ebenen Spannungszustand die Diskussion des Fliessgesetzes in den 
Ecken des aus (1.4) zu gewinnenden Tresca-Sechsecks (Figur 1) zeigt; hier kann 
man indessen diesen Fall, da man die Endlichkeit der iibrigen Verzerrungsge- 
schwindigkeiten postulieren muss, ausschliessen. 

Wenn also eine Achse Hauptachse des Spannungstensors ist, dann ist sie unter 
den in Abschnitt 1 formulierten allgemeinen Voraussetzungen iiber die Fliessbe- 
dingung und das Fliessgesetz auch Hauptachse des Tensors der Verzerrungsge- 
schwindigkeiten, wie das von H. GEIRINGER [7] erstmals festgestellt worden ist. Da- 
gegen ist die Umkehrung dieses Satzes und damit die strenge Richtigkeit der Aus- 
sage, dass die beiden Tensoren koaxial seien, keineswegs evident. Da namlich ein 
symmetrischer Tensor zweiter Ordnung mehr als ein Tripel von Hauptachsen (das 
heisst deren unendlich viele) aufweisen kann, ist es durchaus denkbar, dass der 
Spannungstensor in einem Punkt drei wohldefinierte Hauptachsen, der Tensor der 


A ae 


WN 


Re ee NT 
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Verzerrungsgeschwindigkeiten aber noch weitere Hauptachsensysteme besitzt, so 
dass nicht jede seiner Hauptachsen auch eine solche des Spannungstensors ist. 
Dieser Fall tritt beispielsweise stets dann ein, wenn alle V erzerrungsgeschwindig- 
keiten null sind, das heisst, wenn das zum betreffenden Punkt gehdorende Raum- 


element zwar plastifiziert ist, sich aber wie ein starrer K6rper, bzw. gar nicht bewegt. 


Der Schluss von. (2.2) auf (2.1) wird sich daher, wenn tiberhaupt, nur unter der 
Voraussetzung ziehen lassen, dass sich das Material wirklich plastisch deformiert. 


O 


Figur 1 


Tresca-Sechseck fiir den ebenen Spannungszustand. 


_ Das hier aufgeworfene Problem ist von Bedeutung fiir den ebenen Verschie- 
bungszustand. Dieser wird, wenn 1%, y die Verschiebungsebene ist, durch die Forde- 
rungen 

& = Vy2 = Vex = 0 (2.3) 


samt der Bedingung definiert, dass die iibrigen Verzerrungsgeschwindigkeiten nur 
von * und y abhangen. Aus den letzten beiden Beziehungen (2.3) wird in der 
Literatur allgemein auf (2.1) geschlossen?). Dieser Schluss, so sehr er sich auch auf- 
drangt, folgt aber nicht ohne weiteres aus (2.2). Es ist auch nicht einzusehen, wie 
er sich physikalisch rechtfertigen lasst. Wohl ist man versucht, die Symmetrie 
heranzuziehen und daraus, dass sich zwei benachbarte, zur Ebene ¥, y parallele 
Schichten gleich bewegen, auf die Abwesenheit der die Symmetrie st6renden Schub- 
spannungen t,, und t,, zu schliessen. Symmetriebetrachtungen sind aber in der 
Plastizitatstheorie gefahrlich und entschieden abzulehnen. Bekanntlich liegt es auch 
beim Zugversuch nahe, aus Symmetriegriinden darauf zu schliessen, dass die 
Querkontraktion in allen Querrichtungen gleich gross sei. Aus der Fliessbedingung 
von TreEsca und dem zugehérigen Fliessgesetz folgt aber (Figur 1) nur, dass im 
Punkte A(o,, 0) des FlieBsechsecks fiir den ebenen Spannungszustand der Vek- 
tor é,, é, der Verzerrungsgeschwindigkeit im schraffierten Gebiet liegen muss, und 
da mit Riicksicht auf die Inkompressibilitat des Materials nur noch die Bedingung 
é, + é, + é, = 0 hinzukommt, braucht die Symmetrieforderung é, = é, keineswegs 
erfiillt zu sein. 
Bei Beschrankung auf die Fliessbedingung von v. Mises und das zugehorige 
Fliessgesetz ist der Schluss von (2.2) auf (2.1), sofern man plastischen Fluss voraus- 
setzt, korrekt und leicht zu ziehen (vergleiche zum Beispiel [1], S. 125). Aus (1.3) 


2) Siehe zum Beispiel [7]; [8]; [1], S. 126; [3], S. 268; [9]. 
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erhalt man namlich 0F/0S, = 1 und 0F/0S; = 0, also aus (2.2) nach (1.6) 
Vyg = 2ATy p= 9, Yep = 2ZAtz =O, (2.4) 


und da A = 0 einer starren Bewegung entsprechen wiirde, folgt hieraus (2.1). 

Im folgenden soll dieser Schluss verallgemeinert und nachgewiesen werden, dass 
unter den weitgefassten Voraussetzungen von Abschnitt 1 betreffend Fliessbedin- 
gung und Fliessgesetz beim plastischen Fluss jede Hauptachse des Tensors der Ver- 
zerrungsgeschwindigkeiten auch eine solche des Spannungstensors ist. Im Anschluss 
hieran soll gezeigt werden, dass sich dieses Resultat zwar nicht direkt, aber unter 
einer einschrankenden Voraussetzung auch auf gewisse Kérper mit Verfestigung 
ubertragen asst. 


3. Reduktion des Spannungszustandes 


Bei einer Drehung des Koordinatensystems um die z-Achse (Figur 2) transfor- 


mieren sich die Spannungskomponenten o,,..., T,,,--.. gemass 
o', = 0,cos*y + 2t,, cosysing + a, sin’@, 
o, =o0,sin*g — 27,,, cosp sing + 6, COs* 9; 
Thy = (dy — Gz) COS sin py ==T, , (COS? » — sin? @) , (3.1) 
Th. = Tz, COS P+ T,, Sing, i 
Tye = —Tzy_ SY + Ty, COSY, 
Oo, =6,. 


Die ersten drei Beziehungen (3.1) sind die Transformationen des ebenen Span- 
nungszustandes o,, ¢,, T»,, die sich geometrisch mit dem Spannungskreis interpre- 
tieren lassen; die beiden folgenden sind diejenigen eines Vektors t,,, t,, in der 
Ebene x, y. Entsprechende Transformationen gelten fiir die Komponenten é€,, .. 


<9. 


Yy,/2, ...des Tensors der Verzerrungsgeschwindigkeiten. 
z 
y' 
y 
0 
| Dwmeeare 
x 
Figur 2 


Rotation des Koordinatensystems. 


Wenn die z-Achse Hauptachse des Tensors d indigkei 
ist, dann gelten die beiden Bina (252), wad dwar ate 
auch im gedrehten Koordinatensystem. Mit 4 = 0 oder 0F /OS, = OF. 10s hee 
samtliche Verzerrungsgeschwindigkeiten null, und da man ae Fall ‘der starren 
Bewegung wieder auszuschliessen hat, miissen in (2.2) beide Klammern null sein 
ohne dass gleichzeitig 0F/0S, und 0F /0S, verschwinden. Das ist aber nur méglich, 


ATTA 


a wa Bi ee ee 
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wenn ihre Determinante 

OS rie sell: sa OSes, 9 

O%ye > OF ney 7 ee 3.2) 


zu zx UZ 


ist. Fiihrt man (3.2) mit (1.6) aus, so kommt 


(oy — Oy) Tyz Tee + Tey We — Tee) = 0, (3.3) 


und da die Bedingungen (2.2) auch im gedrehten Koordinatensystem gelten, muss 


(3.3) auch in den gestrichenen Spannungskomponenten erfiillt, mithin unabhangig 
von der Wahl des Winkels 


(0% — Oy) Tye Te + Try (Ty 2 — Teh) = 0 (3.4) 
sein. 

Der dritten Transformation (3.1) ist zu entnehmen, dass man stets einen Winkel 
g so finden kann, dass tj,,, = 0 ist. Da ferner die Fliessbedingung (1.1) von der ersten 
Grundinvariante o,-+ o, + o, des Spannungstensors unabhangig ist, darf man von 
den drei Normalspannungen die gleiche Grésse subtrahieren und daher annehmen, 
dass im gedrehten System von den beiden Normalspannungen oj, , eine beliebige 
null sei. Ist dann die andere ungleich null, so folgt aus (3.4), dass eine der beiden 
Schubspannungen f,,,, t,,, verschwinden, also beispielsweisé Thy = 0, = 0 sein muss. 
Damit ist aber der gegebene auf den ebenen Spannungszustand Oy» Oy Ty, Zuriick- 
gefiihrt. Sind dagegen o;, unf a; gleichzeitig null, so folgt aus den ersten drei Trans- 
formationen, (3.1) bzw. aus dem Spannungskreis, dass auch in jedem anderen Koor- 
dinatensystem mit z als Achse o, = o) = t;,, = 0 ist. Nach der vierten Transforma- 
tion (3.1) kann dann der Drehwinkel gm so gewahlt werden, dass iiberdies 1; ,, = 0 
gilt, und damit ist der Spannungszustand wiederum auf einen ebenen oj, tj, redu- 
ziert. 

Falls also die Bewegung nicht starr und die z-Achse Hauptachse des Tensors der 
Verzerrungsgeschwindigkeiten ist, dann lasst sich der gegebene Spannungszustand 
durch eine Drehung um die z-Achse in jedem Fall auf einen ebenen o,, o,, t,, 
zuriickfiihren. Somit ist 7 eine Hauptachse des Spannungszustandes, und es muss 
jetzt nur noch gezeigt werden, dass auch z eine solche ist. 


4. Der ebene Spannungszustand 


Der Spannungskreis des ebenen Spannungszustandes o,, o,, T,, ist durch die 
Abszisse seines Mittelpunktes und seines Radius, also durch die Gréssen 


I,=o,+ 9,, Fa = (Gy ey)? Aes tp, (4.1) 


vollstandig bestimmt. Daher muss sich die Fliessbedingung (1.1) hier auch in der 
Form 


G(I,, I,) = 0 (4.2) 
und das Fliessgesetz in der Gestalt 
i; = =A [SF + 257 (on — oa]; 
é, = =4[F- -2 37, vo], (4.3) 
Yur =A a = ee 


ZAMP X1I/11 
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anschreiben lassen. In der Tat iiberzeugt man sich leicht davon, dass aus (1.2) mit 


Og = Tze = Vey — 0 


one 


1 
108 


(Fee OT Ta) (4.4) 


il 
<ailld ee 
folgt. Damit ist aber (4.2) gerechtfertigt. 

Nach (4.3) fiihrt die Bedingung 7, = 0, wenn die starre Bewegung wieder aus- 


geschlossen wird, auf ee 
wie == 0), 4.5 
00 ary ne (4.5) 


Da é, und é, frei wahlbar sind, folgt hieraus t,, = 0, denn das Verschwinden des 
zweiten Faktors ware nach (4.3) mit der Beziehung é, = é, gleichbedeutend. Damit 
ist aber nachgewiesen, dass z eine Hauptachse des Spannungszustandes ist. 


5. Kérper mit Verfestiguns 


Bisher wurde angenommen, dass das betrachtete Material idealplastisch sei. 
Lasst man jetzt auch Korper zu, die sich nach der Regel von. PRAGER [10] oder 
nach der vom Verfasser vorgeschlagenen Modifikation ([9], S. 55) derselben ver- 
festigen, so sind in den Beziehungen von Abschnitt 1 die Spannungskomponenten 
rechterhand durchwegs durch die Ausdriicke o,— a,,..., Ty, — Gy,, --- Zu ersetzen, 
und zwar gilt im Falle der Pragerschen Verfestigungsregel beim plastischen 
Fliessen [11] 


: Zeek: ; : 
Og = 3 Chay rer s Shyg = 3 OPyg ares (eu) 
nach ihrer Modifikation dagegen 
Dog = (Ogg — eg) ran 9p Oy, = Py) Ola oasis (5.2) 


Dabei sind ¢ und yw positive Konstanten, deren Bedeutung in diesem Zusammenhang 
unerheblich ist, und ferner gilt zu Beginn des plastischen Flusses 


Ay 60g lo Or. 


US 


Nimmt man an, dass die z-Achse in irgendeinem Stadium des plastischen 
Flusses Hauptachse des Spannungstensors, mithin (2.1) erfiillt sei, so folgt aus den 
in der angegebenen Weise verallgemeinerten Beziehungen (1.5) und (1.6) noch nicht 
(2.2); eine Hauptachse des Spannungstensors ist also im allgemeinen nicht Haupt- 
achse des Tensors der Verzerrungsgeschwindigkeiten. Das hangt damit zusammen, 
dass das Material mit der Verfestigung anisotrop wird. Nimmt man aber die 
Voraussetzung ty, = %» = 0 hinzu, so kommt jetzt (2.2), und zudem schliesst man 
aus (5.1) bzw. (5.2) auf a, = &,, = 0. Folglich ist jede raumfeste Achse, die wahrend 
des ganzen Fliessvorganges eine Hauptachse des Spannungstensors ist, jederzeit 
auch Hauptachse des Tensors der Verzerrungsgeschwindigkeiten und damit auch 
Hauptachse des Verzerrungstensors. 

Nimmt man umgekehrt an, dass die z-Achse in irgendeinem Stadium des pla- 
stischen Flusses Hauptachse des Tensors der Verzerrungsgeschwindigkeiten, mithin 
(2.2) erfiillt sei, so kann man den Spannungszustand o, — « 
gemdss Abschnitt 3 auf einen ebenen Orn hyp (Oh == Coan; ax, a 
im gedrehten Koordinatensystem insbesondere i! 
analog wie in Abschnitt 4, dass DT es 


3 ge Oe 
%,, Zuruckfiihren, wobei 
ze — %z, = 0 wird. Sodann zeigt man 
= 0 sein muss; da aber hieraus und aus 
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Tze — %~ = 9 nicht (2.1) folgt, ist eine Hauptachse des Tensors der Verzerrungs- 
geschwindigkeiten im allgemeinen nicht Hauptachse des Spannungstensors. Nur 
mit der weiteren Voraussetzung «,, = %,, = 0 erhalt man (2.1); da aber dann aus 
(5.1) bzw. (5.2) auch a,, = &,, = 0 folgt, ist jede raumfeste Achse, die wahrend des 
ganzen Fliessvorganges Hauptachse des Tensors der Verzerrungsgeschwindigkeiten 
und somit auch eine solche des Verzerrungstensors ist, jederzeit auch eine Haupt- 
achse des Spannungstensors. 

Damit ist insbesondere auch fiir K6rper, die sich gemass (5.1) oder (5.2) ver- 
festigen, nachgewiesen, dass im Falle des ebenen Verschiebungszustandes (2.3) die 


Normale z zur Verschiebungsebene stets eine Hauptachse des Spannungszustandes 
TStS\ 
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Summary 


In an ideally plastic material any principal axis of the stress tensor is also a 
principal axis of the strain rate tensor. It is shown here that in plastic flow also 
the reverse is true. Besides, the validity of either statement is discussed for mate- 
rials obeying PRAGER’s hardening rule or its modification proposed by the author. 


(Eingegangen: 17. Dezember 1959.) 
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Festschrift Richard Grammel. Zum 70. Geburtstag am 3. Marz 1959. Unter 
Mitwirkung von A.'Berz, K. Krotrer, K. MARGUERRE und E. METTLER heraus- 
gegeben von K. v. SANDEN (Springer-Verlag, Berlin 1959). 372S., 217 Abb.; 
DM 52.-. 

RIcHARD GRAMMEL ist uns vor allem als Griinder und Herausgeber des Inge- 
nieur-Archivs ein Begriff, als Redaktor der Mechanikbande des «alten» H andbuches 
dey Physik, als Autor des bekannten Buches iiber den Kreisel, seine Theorie und seine 
Anwendungen, sowie neben C. B. BrzzENo als Schépfer des berihmten Werkes 
iiber technische Dynamik. 

Die vorliegende Festschrift fiigt sich als Band 28 in den Rahmen des Ingenieur- 
Archivs und enthalt vierzig Aufsatze, die sich iiber das ganze Gebiet der Mechanik 
erstrecken und damit in schénster Weise fiir die Weite des Problemkreises zeugen, 
den GRAMMEL selbst durch seine wissenschaftliche Tatigkeit befruchtet hat. Diese 
Beitrage, von denen die meisten auch den Geist des Lehrers und Freundes Gram- 
MEL dokumentieren, kénnen hier nicht im einzelnen gewiirdigt oder auch nur auf- 
gezahlt werden. Immerhin sei auf das sympathische Geleitwort aus der Feder von 
BIEzENO noch besonders hingewiesen, das mit der folgenden Beurteilung des 
Geehrten schliesst: «Fassen sie [die Schiiler, Freunde und Fachgenossen...] in 
grosster Kiirze zusammen, wobei jedem Wort volle Bedeutung zukommt, so ehren 
sie in ihm den grossen Gelehrten, der seinem Fachgebiet national und international 
sehr bedeutsame Dienste geleistet hat; den vorziiglichen Lehrer; den integren, in 
seinen Auffassungen nicht zu erschiitternden Menschen, der in einer ebenso wahn- 
witzigen wie verbrecherischen Zeit seiner Hochschule, seiner Stadt und seinem Land 
mit Einsatz seiner ganzen Persdnlichkeit und seiner ganzen Arbeitskraft mit 
grossem Erfolg gedient hat.» H. ZIEGLER 


Principles and Applications of Random Noise Theory. Von Jutius 
S. Benpat (John Wiley & Sons, New York and Chapman & Hall, London 1958). 
431 S., 106 Fig.; $11.-. 

Das Bediirtnis, die vielen zerstreuten Erkenntnisse und Methoden iiber die 
Behandlung von gerduschbehafteten Systemen in zweckmassiger Weise zusammen- 
gestellt vor sich zu haben, wird mit dem vorliegenden Werk fast ideal gestillt. Da 
es sich in erster Linie an den Ingenieur wendet, miissen notgedrungen auch die 
statistischen Grundlagen und das mathematische Riistzeug dargestellt werden. So 
behandelt das 1. Kapitel Zufallsprozesse, Charakterisierung der technischen Sy- 
steme und Korrelationsfunktionen. Leistungsspektren, aufbauend auf Fourier- 
Reihe und Fourier-Integral, werden im 2. Kapitel diskutiert und der Zusammen- 
hang mit der Korrelationsfunktion gezeigt, sowohl fiir stationare wie nichtstationare 
Zufallsprozesse. Das 3. Kapitel ist der Wahrscheinlichkeitstheorie gewidmet und 
angewandt auf die Analyse von registrierten Gerduschen. Im 4. Kapitel wird eine 
allgemeine Theorie gegeben iiber die optimale lineare Vorhersage und Filterung 
basierend auf der Methode der kleinsten Quadrate und anwendbar auf lineare 
Systeme mit veranderlichen Koeffizienten, welche durch beliebige nichtstationare 
Zufallsprozesse gestort werden. Die folgenden sechs Kapitel befassen sich mit phy- 
sikalischen und technischen Anwendungen, wobei aber das mathematische Werk- 
zeug entsprechend erganzt wird, Eine besondere Rolle spielt die Autokorrelations- 
funktion, welche eine exponentiell gedampfte Cosinusfunktion ist, zur Darstellung 


CN “nN ay fi 7» ae bi * we \ 
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von vielen vorkommenden Gerauschquellen (5. Kapitel). Fiir solche kann mit Vor- 
teil fiir die Messung von Effektivwerten des Ausgangsfehlers von linearen Systemen 
mit veranderlichen Koeffizienten eine Analogrechenmaschine eingesetzt werden 
(6. Kapitel). Eine wichtige Rolle spielt die Abschatzung der statistischen Fehler bei 
Autokorrelationsmessungen (7. Kapitel). Das 8. Kapitel zeigt, wie Probleme des 
gestorten Radio-Doppel-Empfangs zu behandeln sind, indem hier die Betrachtung 
der Umhiillenden und der Korrelation von Zufallsgerausch dargestellt wird. Dabei 
spielt auch die automatische Verstarkereinstellung eine Rolle, welche beriicksichtigt 
wird. Den optimalen zeitveranderlichen Filtern wird das 9. Kapitel gewidmet unter 
der Voraussetzung, dass das Gerausch durch eine exponentielle Cosinus-Korrela- 
tionsfunktion charakterisiert sei. Das 10. Kapitel befasst sich schliesslich noch mit 
dem schwierigen Problem der Zahl der Nulldurchgange pro Zeitabschnitt, verur- 
sacht durch Gerdusche. Die Betrachtung von Rice wird dahingehend erweitert, 
dass als Signal eine Sinuswelle mit konstanter Frequenz aber statistisch verander- 
licher Amplitude und Phase zugrunde gelegt wird. 

Da das Buch sowohl dem in der Forschung und Entwicklung stehenden Fach- 
mann wie auch dem fortgeschrittenen Studenten dienen will, sind nebst neuesten 
Erkenntnissen auch die grundlegenden Betrachtungen beriicksichtigt worden. Wo 
dies nicht der Fall sein konnte, ist auf die entsprechende Literatur empfehlend hin- 
gewiesen. Ob zu jedem Kapitel eine Einfiihrung und eine Schlussbetrachtung not- 
wendig sei, ist Geschmacksache, sicher ist, dass das vorliegende Buch eine spiirbare 
Liicke in der Handbibhothek zu schliessen vermag fiir alle diejenigen, die mit sta- 
tistisch zufalligen St6rungen in ihren zu bauenden Systemen rechnen miissen. 

H. WEBER 


. Instationare Warmespannungen. Von Hetnz Parkus (Springer-Verlag, 
Wien 1959). 165 S., 34 Abb.; sFr. 38.90. t 

Das Buch ist als Fortsetzung des vom Verfasser zusammen mit Ernst MELAN 
bereits 1953 im gleichen Verlag publizierten Bandchens Wdarmespannungen infolge 
stationdver Temperaturfelder konzipiert. Es geht aber nicht nur in dem Sinne weiter, 
dass jetzt nichtstationére Temperaturfelder betrachtet werden, sondern auch 
darin, dass neben elastischen auch viskoelastische und plastische K6rper in den 
Kreis der behandelten Probleme einbezogen werden. 

Dem Ingenieur, der nach Lésungen in konkreten Fallen sucht, kann das Buch 
bestens empfohlen werden. Es gibt eine auch fiir den Forscher auf diesem oder 
verwandtem Gebiete sehr wertvolle und recht vollstandige Ubersicht tiber den 
heutigen Stand der Forschung und die zur Verfiigung stehenden Lésungen; zudem 
ist die Dokumentation vorbildlich. Als Einfithrung dagegen ist das Buch weniger 
geeignet, da die theoretische Fundierung gerade in den entscheidenden Punkten 
sehr kurz gehalten ist, bzw. fehlt. So enthalt das einleitende Kapitel eine relativ 
breite Behandlung rein mechanischer Prinzipien, die iiberall anderswo nachgelesen 
werden kénnen und im vorliegenden Rahmen von sekundarer Bedeutung sind, 
wihrend bei den in diesem Zusammenhang entscheidenden Grundtatsachen der 
Warmeleitung die Exposition ganz fehlt und durch einige Literaturhinweise er- 
setzt ist. H. ZIEGLER 


Theorie der Stabilitat einer Bewegung. Von Jorr G. Marky, in deutscher 
Sprache herausgegeben von WoLFrGaNnG HAHN und Rotr Retssic (R. Oldenbourg, 
Miinchen 1959). 402 S., 20 Abb.; DM 47.-. 

Die Entwicklung der Regelungstechnik, der Astronautik und vieler anderer 
Gebiete bringt es mit sich, dass der Theorie der Stabilitat von Bewegungen immer 


166 Buchbesprechungen — Book Reviews — Notices bibliographiques ZAMP 
gréssere Bedeutung zukommt. Die klassischen Methoden, die auch heute noch eine 
grosse Rolle spielen, stiitzen sich auf die Linearisierung der Bewegungsditferential- 
gleichungen. Es ist indessen schon friih erkannt worden, dass dieses Verfahren 
logisch nicht einwandfrei ist und in gewissen Fallen ein falsches Bild von den 
wirklichen Stabilitatsverhdltnissen gibt. a 

Die Grundlage fiir eine mathematisch strenge Behandlung des Stabilitats- 
problems ist 1892 in seiner Dissertation durch Lyapunow geschaffen worden. Auf 
dieser Basis hat sich die seitherige Forschung entwickelt; dabei hat vor allem die 
russische Schule grosse Erfolge erzielt. Das vorliegende Werk ist von einem kompe- 
tenten Kenner und Foérderer der Materie in streng mathematischer Form ge- 
schrieben und fasst die vielen Einzelarbeiten auf diesem Gebiet zu einem ausge- 
zeichneten Lehrbuch zusammen. Die russische Fassung datiert von 1952; der 


deutsche Text enthalt aber viele Verbesserungen von seiten des Verfassers wie der 
Ubersetzer. H. ZIEGLER 


Gasentladungsréhren in der Nachrichtentechnik. Nachrichtentech- 
nische Fachberichte NFT, Bd. 9. Herausgegeben von J. WosNiIkK (Friedr. 
Vieweg & Sohn, Braunschweig 1957). 62 S., 99 Fig.; DM 8.50. 

Das vorliegende Heft enthalt zwolf Vortrage, die von europdischen Gasentla- 
dungs-Spezialisten an der Diskussionstagung in Aachen (Februar 1957) gehalten 
worden sind. Sie sind darum interessant, weil sie in einen Zeitpunkt fallen, in dem 
sich das allgemeine Interesse der Elektroniker von den Gasentladungsro6hren eher 
abwandte, um sich mehr auf die Halbleitertechnik zu konzentrieren. Die Grinde 
dieses Verhaltens — vor allem der Automatiker — liegen zum Teil in den fiir elek- 
tronische Schaltelemente relativ langen Schaltzeiten, in Unannehmlichkeiten im 
Zusammenhang, mit dem Léschvorgang sowie in den hohen Betriebs- und Steuer- 


spannungen. Die vorteilhaften Eigenschaften dieser Bauelemente, wie hohe Le- — 


bensdauer und grosse Robustheit im Betrieb, rechtfertigen offenbar weitere Ent- 
wicklungsarbeiten zur Verbesserung der Rohreneigenschaften. 

So wird berichtet tiber Massnahmen zur Verkiirzung und Stabilisierung der 
Schaltzeit, tiber Untersuchungen von Ziindempfindlichkeit und Léschzeit. Im 
Blick auf tibertragungstechnische Anwendungen interessieren die Studien tiber 
Stérschwingungen und deren Vermeidung, wahrend der Schaltungsingenieur niitz- 
liche Hinweise fiir die Entwicklung von Tor- und Zahlschaltungen findet. Ein 
besonderer Abschnitt ist den Sperréhren gewidmet, die im Frequenzbereich bis 
10000 MHz als Schaltdioden arbeiten miissen und daher besondere physiko- 
chemische Probleme aufwerfen. Der letzte Beitrag behandelt die Glimmroéhren als 
Uberspannungsableiter zum Schutz von Fernmeldeanlagen. 

Alle Abschnitte sind durch gute graphische Darstellungen, Oszillogramme, 
Konstruktionsskizzen und Literaturangaben bereichert. Das Heft bietet einen 
guten Uberblick tiber den Stand der Entwicklung und die Moglichkeiten der tech- 
nischen Anwendung der Gasentladungsréhren und diirfte Physiker und Ingenieure 
in gleicher Weise interessieren. F.KuMMER 


Praktische Spannungsoptik. Von L.Féppr und E. Moncu (Springer- 
Verlag, Berlin 1959). 209 S., 163 Abb.; DM 30.-., 

Die zweite Auflage dieses Buches, das in erster Linie fiir den in der Praxis tati- 
gen Ingenieur geschrieben ist, wurde zum grossen Teil vollstandig neu bearbeitet, 


da seit dem Erscheinen der ersten Auflage verschiedene Methoden der Spannungs- 
optik erheblich verbessert worden sind. 
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Der erste Abschnitt gibt die theoretischen und experimentellen Grundlagen 
der ebenen und raumlichen Photoelastizitat sowie eine Behandlung der mechani- 
schen Ahnlichkeitsgesetze, die fiir die Ubertragung der Ergebnisse vom Modell 
auf die Wirklichkeit wichtig~sind. Im zweiten Abschnitt sind die besonderen 
Verfahren kurz beschrieben, wahrend der dritte einige typische Anwendungen 
enthalt, die am Spannungsoptischen Laboratorium der Technischen Hochschule 
Miinchen ausgefiihrt worden sind. 

Im Gegensatz zur ersten Auflage, in welcher aus verschiedenen Griinden die 
einfache spannungsoptische Apparatur (mit Filter, ohne Linsen und mit diffusem 
Licht) der tiblichen optischen Bank bevorzugt wurde, werden jetzt beide Ein- 
richtungen einander gegeniibergestellt. Das neueste in Miinchen konstruierte 
Gerat vereinigt die Vorteile der ersteren (leichte Handhabung, gute Beobachtungs- 
méglichkeit grosser Flachen des Modells) mit den Vorteilen der letzteren (paralleler 
Strahlengang, Méglichkeit des Projizierens). 

Das Buch ist ein Zeugnis der jahrzehntelangen Erfahrung beider Verfasser 
auf dem Gebiet der Photoelastizitat und kann daher allen empfohlen werden, 
die sich mit dieser Wissenschaft und deren Anwendung beschaftigen wollen. 

W. SCHUMANN 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLUGcGe und E. Creutz, Band 44, 1.Teil: Instyvuwmentelle Hilfsmittel des Kern- 
phystkers (Springer-Verlag, Berlin 1959). 473 S., 225 Fig.; DM 125.-. 

Der vorliegende Band der Fliggeschen Handbuchreihe befasst sich mit 
Apparaten zur Beschleunigung geladener Teilchen. Das Gebiet des Acceleratoren- 
baues hat sich im Laufe der letzte Jahre zu industriellen Ausmassen ausgeweitet 
und hegt heute nicht mehr ausschliesslich in den Handen des Kernphysikers, 
sondern erdffnet auch Ingenieuren, Mathematikern und Technikern ein weites 
Tatigkeitsfeld. Neuere Arbeitsprinzipien, wie etwa die Entdeckung des Prinzips 
der starken Fokussierung, hatten zur Folge, dass immer hoéhere Werte der Maxi- 
malenergie angestrebt werden, was natiirlich mit einem entsprechenden Anstieg 
des Aufwandes verbunden ist. Es ist daher charakteristisch, dass diejenigen 
Physiker, die gréssere Acceleratoren fiir kernphysikalische Experimente benutzen, 
im allgemeinen nicht mehr in der Lage sind, die Finessen ihres Beschleunigers im 
Detail zu iiberblicken. Und gerade fiir diese Physiker diirfte das vorliegende Buch 
von Nutzen sein: es vermittelt Kenntnisse tiber die wichtigsten Acceleratoren- 
typen und deren Wirkungsweise, ohne dabei auf alle diejenigen Einzelheiten ein- 
zugehen, die fiir den Konstrukteur von Bedeutung sind. 

Die zwei ersten Kapitel befassen sich mit dem Kaskadengenerator (E.BaLpIN- 
GER) und dem Van-De-Graff-Generator (R.G.HERB) und erganzen sich insofern 
vorteilhaft, als die Behandlung des Kaskadengenerators sich auf die Erzeugung 
der Hochspannung beschrankt, wahrend Probleme, die beiden Acceleratoren 
gemeinsam sind, wie Ionenquellen, Strahlrohre usw., im zweiten Kapitel disku- 
tiert werden. Die Zirkularbeschleuniger werden in vier weiteren Kapiteln behan- 
delt: Zyklotron und Synchrozyklotron (B.L.CouEN), Elektron-Synchrotron 
(R.R.Witson), Betatron (D.W.Kerst) nnd Proton-Synchrotron (G.K.GREEN). 
Die Probleme der Teilchenbahnen, die Bahn- und Phasenstabilitat, der Injek- 
tions- und Extraktionsmechanismus beanspruchen ihrer Bedeutung gemass hier 
den weitesten Raum. Daneben werden aber auch die wichtigsten konstruktiven 
und technischen Fragen nicht vernachlassigt. Wie an allen Stellen des Buches 
werden zur Illustration in der Praxis ausgefiihrte Maschinen besprochen. Ein 
weiteres von SmitH verfasstes Kapitel iiber lineare Acceleratoren wird mit einer 
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allgemeinen Diskussion der Bewegungsgleichung eingeleitet. Elektronen- und 
Ionenbeschleuniger sind dabei separat behandelt. Besonders sorgfaltig werden 
die fiir den Bau so wichtigen Toleranzforderungen besprochen und als Beispiele 
der Elektronenbeschleuniger in Stanford und der Beschleuniger fur schwere Jonen 
in Berkeley angefiihrt. Im letzten Kapitel (D.J.HuGHEs) wird der Reaktor als 
Neutronenquelle fiir kernphysikalische Experimente besprochen. Die Wahl der 
in diesem Kapitel diskutierten Objekte erscheint auf den ersten Blick etwas 
willkiirlich. Ein Vergleich mit vorangehenden Banden dieser Handbuchreihe 
zeigt jedoch, dass hier eine Fiille von wertvollen Erganzungen zu anderen Artikeln 
(vgl. zum Beispiel Instrwmentelle Hilfsmittel des Kernphystkers, 2.Teil) geboten 
wird. Man findet Angaben iiber Energiespektren von Reaktorneutronen, charak- 
teristische Neutronenmessanordnungen sowie einen kurzen Uberblick tiber Strah- 
jungsschutzmassnahmen. 

Wenn man den Zweck des vorliegenden Buches in der eingangs skizzierten 
Weise versteht und nicht eine erschépfende Darstellung der behandelten Gegen- 
stande erwartet, wird man —sofern nicht der hohe Preis davon abschreckt — gern 
zum Erwerb dieses Bandes schreiten. F. HEINRICH 


An Introduction to Plasticity. Yon WiLL1am PRaAGER (Addison-Wesley 
Publ. Co., Inc., Reading 1959). 148 S., 82 Fig.; $6.50. ; 

Das Buch ist als Ubersetzung des Bandchens Probleme der Plastizitdtstheorie 
angekiindigt, das 1955 im Verlag Birkhauser erschienen und in ZAMP 7, 471 (1956) 
besprochen worden ist. Tatsachlich stellt es eine modernisierte und erweiterte 
Fassung der deutschen Ausgabe dar. Neu sind die Abschnitte tiber Eindeutigkeit 
des anfanglichen Flusses im starrplastischen Medium unter gegebenen Lasten, tiber 
optimale Formgebung (Limit Design) und iiber einen Fall von modifiziert-pseudo- 
stationdrem Fluss. Auch in den tibrigen Abschnitten (zum Beispiel bei den Platten 
und in der Behandlung des Einflusses von Anderungen in der Geometrie des sich 
_ verformenden KGrpers) sind neuere Resultate beriicksichtigt. Die Literaturhinweise 
sind stark vermehrt sowie dem heutigen Stand angepasst worden, und schliesslich 
enthalt die englische Fassung eine gréssere Zahl von dusserst interessanten Ubungs- 
aufgaben. 

Man kann sich fragen, ob der neue Titel dem Inhalt des Buches gerecht wird. 
Einfithrungen in die Plastizitat gibt es viele; als knappe, meisterhaft klare, an- 
schauliche und zudem ausgezeichnet dokumentierte Ubersicht iiber die zurzeit 
wichtigsten Probleme und Verfahren steht dieses Werk aber allein. H. ZIEGLER 


Approximate Methods of Higher Analysis. Von L. V. Kanrorovicu und 
V. I. KryLov (P. Noordhoff Ltd., Groningen 1958). 681 S., 68 Fig.; $ 17.-. 

Der vorliegenden englischen Ubersetzung liegt eine vierte, offenbar unveran- 
derte Auflage des russischen Originals zugrunde. Die Autoren bemerken prinzipiell, 
dass heute ein viel umfangreicheres Erfahrungsmaterial iiber die Anwendungen der 
behandelten Methoden zur Verfiigung steht als zur Zeit des ersten Erscheinens 
dieses Buches und dass der Einsatz von Rechenautomaten zum Teil neue Gesichts- 


punkte gebracht hat, und sie behalten sich fiir die Zukunft eine entsprechende 
_ totale Neubearbeitung des Stoffes vor. : 


Der Inhalt umfasst: Lésung von partiellen Differentialgleichungen durch 


Reihenansatze ; approximative Lésung von Fredholmschen Integralgleichungen; 
Differenzenmethoden ; Methoden der Variationsrechnung; Konforme Abbildungen, 
mit Anwendung auf harmonische und biharmonische Probleme; die Schwarzsche 


Methode (Dirichlet- und Neumann-Problem), P. LAucHLI 
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Neuerscheinung. 


EXPERIENTIA SUPPLEMENTUM VIII 


Bericht tiber die Tatigheit der Schweizerischen 
Studienkommission fiir Atomenergie von 1946 bis 1958 


112 Seiten mit 45 Abbildungen. Broschiert Fr./DM 9.— (1960) 


INHALT 
Geschichte und Tatigkeit der Schweiz. Studienkommission fiir Atomenergie 1946-1958 von Dr. A. Krethlow, Bern. 


BASEL: Physikalisches Institut der Universitat Basel, von Prof. Dr. P. Huber, Physikalisch-Chemisches Institut der 
Universitat Basel, von Prof. Dr. W. Kuhn. 


BERN: Chemisches Institut der Universitat Bern, von Prof. Dr. W. Feitknecht und Prof. Dr. W. Buser. Medizinisch- 
Chemisches Institut der Universitat Bern, von Prof. Dr. med. H. Aebi. Physikalisches Institut der Universitat Bern, 
von Prof. Dr. F. G. Houtermans. Theodor-Kocher-Institut der Universitat Bern, von Prof. Dr. J. Mattauch und — 
Dr. H. Hintenberger. Réntgeninstitut des Inselspitals Bern, von Prof. Dr. A. Zuppinger. Arbeitsausschuss fiir die 
Untersuchung schweizerischer Mineralien und Gesteine auf Atombrennstoffe und seltene Elemente, von Prof. Dr. . 
F. de Quervain und Prof. Dr. Th. Hiigt. 


GENF: Physikalisches Institut der Universitat Genf, von Prof.Dr. R. Extermann. Physikalisches Institut der Univer- 
sitat Genf, von Prof. Dr. E. Stueckelberg. 


LAUSANNE: Laboratoire de Recherches Nucléaires, Ecole Polytechnique de l’Université de Lausanne, von Prof. 
Ch. Haenny. 


NEUENBURG:; Physikalisches Institut der Universitat Neuenburg, von Prof. Dr. A. iia! und Dr. R. Payot. 
Physikalisches Institut der Universitat Neuenburg, von Prof. Dr. J. Rossel. S 


ZURICH: Physikalisches Institut der Eidgendssischen Technischen Hochschule, von Prof. Dr. P. Scherrer. ceiaages 
rium fiir Anorganische Chemie der Eidgendssischen Technischen Hochschule, Ziirich, von Prof. Dr. W. Treadwell. 
Chemisches Institut der Universitat Ziirich, von Prof. Dr. P. Karrer. Physikalisch-Chemisches Institut der Universitat 
Zarich, von Prof. F © K. Clusius. Physikalisches Institut der Universitat Ztirich, von Prof. Dr. Hans H. Staub. 


In der vorliegenden Publikation wird in gedrangter Form iiber die Tatigkeit der Schweizerischen 
Studienkommission fiir Atomenergie (SKA) wahrend der Jahre 1946 bis 1958 berichtet, sowie tiber | 
die Forschungsarbeiten der mit ihr zusammenarbeitenden Hochschulinstitute. 

Die Studienkommission fiir Atomenergie wurde aus der Erkenntnis heraus, dass die Schweiz hin- 
sichtlich der Forschungen auf dem Gebiete der Kernphysik nicht abseits stehen kann und darf, auf 
Grund eines Bundesbeschlusses vom Jahre 1946 geschaffen. Sie wurde zum planenden und ausfiih- 
renden Organ bestimmt fiir die vom Bund tibernommene Aufgabe der Forderung der Forschung auf 
dem Gebiet der Atomenergie, wobei der Bund die notwendigen finanziellen Mittel zur Verfiigung 
stellte. Hiermit war es méglich, die Forschungsarbeiten in der Schweiz auf eine breitere Basis zu 
stellen und auch in vermehrtem Masse Kernphysiker und Spezialisten auf den angewandten Zweigen 
der Kernphysik heranzubilden. Samtliche schweizerischen Physikalischen sowie auch einige Physi- 
kalisch-Chemische und Medizinische Institute erklarten sich zur Mitarbeit bereit. Uber die von ihnen 
erhaltenen interessanten Forschungsergebnisse wird in den einzelnen Abschnitten referiert. 

Die Publikation stellt eine wertvolle Zusammenfassung der an unseren Universitaten und der Eidg. 
Techn. Hochschule geleisteten Forschungstatigkeit auf dem Gebiete der Atomenergie wahrend der 


Jahre 1946 bis 1958 dar. 
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